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1 Introduction 

Constructing solutions for (systems of) partial differential equations is obviously difficult - in 
particular for non-linear systems. Elie Cartan ^ proposed to construct first infinitesimal solu- 
tions or integral elements. These are possible tangent spaces to (prolonged) solutions. Thus they 
always lead to a linearisation of the problem and their explicit construction requires essentially 
only straightforward linear algebra. In the Cartan-Kahler theory [3l|T2], differential equations 
are represented by exterior differential systems and integral elements consist of tangent vectors 
pointwise annihilated by differential forms. 

Vessiot j28j proposed in the 1920s a dual approach which does not require the use of exterior 
differential systems. Instead of individual integral elements, it always considers distributions 
of them generated by vector fields and their Lie brackets replace the exterior derivatives of 
differential forms. This approach takes an intermediate position between the formal theory of 
differential equations [13 \ \19 \ [23] and the Cartan-Kahler theory of exterior differential systems. 
Thus it allows for the transfer of many techniques from the latter to the former one, although 
this point will not be studied here. 

Vessiot's approach may be considered a generalisation of the Frobenius theorem. Indeed, if 
one applies his theory to a differential equation of finite type, then one obtains an involutive 
distribution such that its integral manifolds are in a one-to-one correspondence with the smooth 
solutions of the equation. For more general equations, Vessiot proposed to "cover" the equation 
with infinitely many involutive distributions such that any smooth solution corresponds to an 
integral manifold of at least one of them. 

Vessiot's theory has not attracted much attention: presentations in a more modern language 
are contained in [Sj [24]; applications have mainly appeared in the context of the Darboux 
method for solving hyperbolic equations, see for example [27]. While a number of textbooks 
provide a very rigorous analysis of the Cartan-Kahler theory, the above mentioned references 
(including Vessiot's original work [28]) are somewhat lacking in this respect. In particular, the 
question of under what assumptions Vessiot's construction succeeds has been ignored. 

The purpose of the present article is to close this gap and at the same time to relate the Vessiot 
theory with the key concepts of the formal theory like formal integrability and involution (we 
will not develop it as a dual form of the Cartan-Kahler theory, but from scratch within the 
formal theory). We will show that Vessiot's construction succeeds if, and only if, it is applied to 
an involutive system of differential equations. This result is of course not surprising, given the 
well-known fact that the formal theory and the Cartan-Kahler theory are equivalent. However, 
to our knowledge an explicit proof has never been given. As a by-product, we will provide 
a new characterisation of integral elements based on the contact map, making also the relations 
between the formal theory and the Cartan-Kahler theory more transparent. Furthermore, we 
simplify the construction of the integral distributions. Up to now, quadratic equations had to 
be considered (under some assumptions, their solution can be obtained via a sequence of linear 
systems. We will show how the natural geometry of the jet bundle hierarchy can be exploited 
for always obtaining a linear system of equations. 

This article contains the main results of the first author's doctoral thesis [6] (a short summary 
of the results has already appeared in [?])• It is organised as follows. The next three sections 
recall the needed elements from the formal theory of differential equations and provide our new 
characterisation of transversal integral elements. The following two sections introduce the key 
concepts of Vessiot's approach: the Vessiot distribution, integral distributions and flat Vessiot 
connections. Two further sections discuss existence theorems for the latter two based on a step- 
by-step approach already proposed by Vessiot. As the proofs of some results are fairly technical, 
the main text contains only an outline of the underlying ideas and full details are given in three 
appendices. The Einstein convention is sometimes used to indicate summation. 
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2 The contact structure 

Before we outline the formal theory of partial differential equations, we briefly review its under- 
lying geometry: the jet bundle and its contact structure. Many different ways exist to introduce 
these geometric constructions, see for example I18[ [2U]. Furthermore, they are discussed in 
any book on the formal theory (see the references in the next section). 

Let TT : £ ^ X he a smooth fibred manifold. We call coordinates x = (x* : 1 < i < n) of 
independent variables and fibre coordinates u = (n" ■ 1 < a < m) in £ dependent variables. 
Sections a : X ^ £ correspond locally to functions u = s(x). We will use throughout a "global" 
notation in order to avoid the introduction of many local neighbourhoods even though we mostly 
consider local sections. 

Derivatives are written in the form = d^^^u" /dx^^ ■ ■ ■ dxn" where fi = (/ii, . . . , is 
a multi-index. The set of derivatives u'^ up to order q is denoted by u^''^; it defines a local 
coordinate system for the g'-th order jet bundle Jgvr, which may be regarded as the space of 
truncated Taylor expansions of functions s. 

The hierarchy of jet bundles Jqir with q = 0,1,2,... possesses many natural fibrations 
which correspond to "forgetting" higher-order derivatives. For us particularly important are 
'^q-i ■ '^i'^ ~^ Jq-i'^ Siud tt'^ : JqTT X. To cach section a : X ^ £^ locally defined by 
(t(x) = (x, s(x)) , we may associate its prolongation jqa : X Jgvr, a section of the fibration n'' 
locally given by jqa(x) = (x, s(x), 9xs(x), 9xxs(x), . . . ) . 

The geometry of the jet bundle Jqir is to a large extent determined by its contact structure. It 
can be introduced in various ways. For our purposes, three different approaches are convenient. 
First, we adopt the contact codistribution C T*{Jq7r); it consists of all one- forms such that 
their pull-back by a prolonged section vanishes. Locally, it is spanned by the contact forms 

n 

= dn^ - ^ < 1^1 < gr, 1 < a < m. 

i=l 

Dually, we may consider the contact distribution Cq C T(Jq7r) consisting of all vector fields 
annihilated by . A straightforward calculation shows that it is generated by the contact fields 

m 

a=10<|Aj|<g 

C^ = a„2, 1^1 =g, l<a<m. (1) 

Note that the latter fields, C^, span the vertical bundle Vtt'^_i of the fibration iTg_i. Thus 
the contact distribution can be split into Cq = Vir'^^i TC. Here the complement Ti is an n- 
dimensional transversal subbundle of T{JqTT) and obviously not uniquely determined (though 
any local coordinate chart induces via the span of the vectors C^''^ one possible choice). Any 
such complement Ti may be considered the horizontal bundle of a connection on the fibred 
manifold vr*^ : JqTT — > X (but not for the fibration TTq_i). Following Fackerell [5], we call any 
connection on tt^ the horizontal bundle of which consists of contact fields a Vessiot connection 
(in the literature the terminology Cartan connection is also common, see for example |15j). 

For later use, we note the structure equations of the contact distribution. The only non- 
vanishing Lie brackets of the vector fields ([1]) are 

[Cr^NCi")] =a.., \v\=q-l. (2) 



Note that this observation implies that the vertical bundle Vtt'^_'^ is involutive. 
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As a third approach to the contact structure we consider, following Modugno [T7j, the contact 
map (see also [SI [23]). It is the unique map Tg : Jgir x TX T(Jg„i7r) such that the diagram 



X 

J vr X TX ^ TUq-lTT) 

X ^ 

((i<jO")oTA-)><^dTA\ ^/TXiq-lo") 

TX 

commutes for any section a. Because of its linearity over vr^^^, we may also consider it a map 
Tq : Jgvr — > T*X ^ T(Jq„i7r) with the local coordinate form 

Jq-ITT 

r, : (x,u(«)) ^ (x,u(«-i);dx* {d,^ + W^^+iA^,)) ■ (3) 

Obviously, Tq{p,d^i) = rp7r^_^(cf^^) and hence {Cg)p = (Tpvr^.^^)-^ (imr,(/9)) for any point 
p € JgTT. Note that all vectors in the image of Tq[p) are transversal to the fibration vr^Zj- 

One of the main applications of the contact structure is given by the following proposition 
(for a proof, see [Gj Proposition 2.1.6] or |23^ Proposition 2.2.7]). It characterises those sections 
of the jet bundle tt'^ : J^vr — > X which are prolongations of sections of the underlying fibred 
manifold it: £ ^ X. 

Proposition 1. A section 'j : X ^ Jgir is of the form 7 = jqa for a section a : X ^ £ if, and 
only if, imrg(7(x)) = T^(^x)''^'!i-i{T^{x) im7) for all points x G X where 7 is defined. 

Thus for any section a : X ^ £ the equality imrq_|_i(jg+i(T(x)) = im.Tx{jqCr) holds and we 
may say that knowing the {q + l)-jet jq^ia{x) of a section a at some x € X is equivalent to 
knowing its g-jet p = jq(T{x) at x as well as the tangent space Tpiivajqu) at this point. This 
observation will later be the key for the Vessiot theory. 

3 The formal theory of differential equations 

We are now going to outline the formal theory of partial differential equations to introduce the 
basic notation. Our presentation follows [23]; other general references are [13 ^ 114 ^ [19]. 

Definition 1. A differential equation of order g is a fibred submanifold TZq C J^tt locally 
described as the zero set of some smooth functions on Jgvr: 

. r <^^(x,uW) = o, 

I (r=l,...,t). 

Note that we do not distinguish between scalar equations and systems. 

We denote by i : TZq ^ Jqir the canonical inclusion map. Differentiating every equation 
in the local representation ^ leads to the prolonged equation TZq+i C Jg+i7r defined by the 
equations (^"^ = and Di^'^ = where the formal derivative Di is given by 

A<^^(x,u(^-^)) = ^(x,u(^))+ E^(-'"^^V.V. (5) 

o<|mI<ij"=i ^ 

Iteration of this process gives the higher prolongations TZq+r ^ Jq+r"^- A subsequent projec- 
tion leads to Tzi^^ = ■Kq'^^ (TZq-^i) C TZq^ which is a proper submanifold whenever integrability 
conditions appear. 
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Definition 2. A differential equation TZq is formally integrable if at any prolongation order 
r > the equality 7^^^^ = TZg+r holds. 

In local coordinates, the following definition coincides with the usual notion of a solution. 
Definition 3. A solution is a section a : X ^ £ such that its prolongation satisfies imjqa C TZ^. 

For formally integrable equations it is straightforward to construct order by order formal 
power series solutions. Otherwise it is hard to find solutions. A constitutive insight of Cartan 
was to introduce infinitesimal solutions or integral elements at a point p € TZq as subspaces 

Q TpTZq which are potentially part of the tangent space of a prolonged solution. 

Definition 4. Let TZq C J^vr be a differential equation and l : TZq Jqir the canonical inclusion 
map. Let IlTZq] = {i*Cq)diH be the differential ideal generated by the pull-back of the contact 
codistribution on TZq (algebraically, I[TZq] is then spanned by a basis of i*Cq and the exterior 
derivatives of the forms in this basis). A linear subspace Up C TpTZq is an integral element at 
the point p € TZq, if all forms in {T\TZq\)p vanish on it. 

The following result provides an alternative characterisation of transversal integral elements 
via the contact map. It requires that the projection t^\'^^ : T^g+i — > T^q is surjective. 

Proposition 2. Let TZq he a differential equation such that TZ^^ = TZq. A linear subspace 
Up C TpTZq such that TpL{Up) lies transversal to the fibration vr^_i is an integral element at the 
point p G TZq if, and only if, a point p G 7^q+i exists on the prolonged equation 7^q+i such that 
^g'^^P) = P and Tpi{Up) C {Ta.Tq+i{p). 

Proof. Assume first that Up satisfies the given conditions. It follows immediately from the 
coordinate form of the contact map that then firstly Tpi{Up) is transversal to 7rg_]^ and secondly 
that every one-form oj G i*Cq vanishes on Up, as imVqj^i^p) C {Cq)p. Thus there only remains 
to show that the same is true for the two-forms du) G L*{dCq ). 

Choose a section 7 : TZq — > 7^g+i such that j{p) = p and define a distribution V of rank 
n on TZq by setting Ti{T>p) = \mTqj^i[^{p)^ for any point p G TZq. Obviously, by construction 
Up C T)p. It follows now from the coordinate form ([3]) of the contact map that locally the 
distribution V is spanned by n vector fields Xi such that = c\'^^ + jp^i.Ca where the 
coefficients 7" are the highest-order components of the section 7. Thus the commutator of two 
such vector fields satisfies 

The commutators on the right hand side vanish whenever pi = or pj = 0, respectively. 
Otherwise we obtain —du" and — , respectively. But this fact implies that the two sums 

fj, 1^ M 

on the right hand side cancel each other and we find that G Cq. Thus we find for 

any contact form w G Cg that 

i*{duj){Xi,Xj) = d(jj{L^,Xi,L^Xj) = i^Xi[u!{L^:Xj)) — i*Xj(a;(i*Xj)) + a;(/-*([Xj, Xj])) . 

Each summand in the last expression vanishes, as all appearing fields are contact fields. Hence 
any form u G L*{dCq) vanishes on V and in particular on Up QT)p. 

For the converse, note that any transversal integral element Up C TpTZq is spanned by linear 
combinations of vectors Vi such that Tpi{vi) = c\'^\p + JpiCtilp where 7^^ are real coefficients. 
Now consider a contact form cj" with = q — 1. Then = dx^ A du'^^i.. Evaluating 
the condition i*{duj'^)\p{vi,Vj) = duj{TpL{vi),Tpi{vj)) = yields the equation 7^?+i;j = 7^?+i^.^j- 
Hence the coefficients are of the form 7^^^^ = 7^_|_i. and a section a exists such that p = jq(j{x) 
and Tpimijqa) is spanned by the vectors Tpi{vi), . . . ,Tpi{vn). This observation implies that Up 
satisfies the given conditions. ■ 
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For many purposes the purely geometric notion of formal integrability is not sufficient, and 
one needs the stronger algebraic concept of involution. This concerns in particular the derivation 
of uniqueness results but also the numerical integration of overdetermined systems |21] . An 
intrinsic definition of involution is possible using the Spencer cohomology (see for example [22] 
and references therein for a discussion). We apply here a simpler approach requiring that one 
works in "good", more precisely: 5-regular, coordinates x. This assumption represents a mild 
restriction, as generic coordinates are 5-regular and it is possible to construct systematically 
"good" coordinates - see [TT]. Furthermore, it will turn out that the use of 5-regular coordinates 
is essential for Vessiot's approach. 

Definition 5. The (geometric) symbol of a differential equation TZq is Afq = Vtt^_^\ti^ D TlZq. 

Thus, the symbol is the vertical part of the tangent space to TZq. Locally, Nq consists of all 
vertical vector fields Yl^=iYli\^\=q''^'^'^u<^ where the coefficients satisfy the following linear 
system of algebraic equations: 

= 0. (6) 

The matrix of this system is called the symbol matrix Mq. The prolonged symbols Afq+r are the 
symbols of the prolonged equations TZq+r with corresponding symbol matrices Mq+r- 

The class of a multi-index ^ = (fii, . . . ,fj,n), denoted els//, is the smallest k for which fi^ is 
different from zero. The columns of the symbol matrix ([6]) are labelled by the v'^^. We order them 
as follows. Let a and P denote indices for the dependent coordinates, and let /i and ly denote 
multi-indices for marking derivatives. Derivatives of higher order are greater than derivatives 
of lower order: if \fi\ < \u\, then -< Uy. If derivatives have the same order = then we 

distinguish two cases: if the leftmost non-vanishing entry m. ^ — v \s positive, then -< Up\ and 

if fi = u and a < P, then -< u^- This is a class-respecting order: if = |z^| and els fi < els u, 

then u'^ -< Uu- Any set of objects indexed with pairs (q, fx) can be ordered in an analogous way. 
This order of the multi-indices /i and v is called the degree reverse lexicographic ranking, and 
we generalise it in such a way that it places more weight on the multi-indices fi and z/ than on 
the numbers a and /? of the dependent variables. This is called the term- over-position lift of the 
degree reverse lexicographic ranking. 

Now the columns within the symbol matrix are ordered descendingly according to the degree 
reverse lexicographic ranking for the multi-indices /.i of the variables in equation ([6]) and 
labelled by the pairs {a,fi). (It follows that, if v'^ and are such that cls/i > clsz^, then the 
column corresponding to v'^ is left of the column corresponding to v^.) The rows are ordered in 
the same way with regard to the pairs (a, fi) of the variables u'^ which define the classes of the 
equations «P'^(x, u^'^)) = 0. If two rows are labelled by the same pair (a,/i), it does not matter 
which one comes first. 

We compute now a row echelon form of the symbol matrix. We denote the number of rows 
where the pivot is of class k by I3q''\ the indices of the symbol Mq, and associate with each 
such row its multiplicative variables x^, . . . ,x^. Prolonging each equation only with respect to 
its multiplicative variables yields independent equations of order q -\- 1, as each has a different 
leading term. 

Definition 6. If prolongation with respect to the non-multiplicative variables does not lead to 
additional independent equations of order g + 1, in other words if 

n 

rankM,+i = (7) 
k=l 
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then the symbol Nq is involutive. The differential equation TZq is called involutive, if it is formally 
integrable and its symbol is involutive. 

The criterion ([7]) is also known as Cartan's test, as it is analogous to a similar test in the 
Cartan-Kahler theory of exterior differential systems. We stress again that it is valid only in 
5-regular coordinates (in fact, in other coordinate systems it will always fail). 



4 The Cartan normal form 

For notational simplicity, we will consider in our subsequent analysis almost exclusively first- 
order equations TZi C J^tt. At least from a theoretical standpoint, this is not a restriction, as any 
higher-order differential equation TZq can be transformed into an equivalent first-order one (see 
for example |23t Appendix A.3]). For these we now introduce a convenient local representation. 



Definition 7. For a first-order differential equation TZi the following local representation, a spe- 
cial kind of solved form. 



u. 



,° = C(x,u^u;,4) 



" i=C-i(x,u^uJ,u^_i) 



1 < a < 
I <j <n, 

/3j"^ <6<m, 

1 < j < n - 1, 

/3j"~^^ <6<m, 



(8b) 



< = </.?(x,n^nf) (8c) 

is called its Cartan normal form. The equations of zeroth order, = (j)'^{^,u^), are called 
algebraic. The functions (p'^ are called the right sides of TZi. (If, for some 1 < k < n, the number 

(k) (k) A 

of equations is (31 = m, then the condition /3J < 5 < m is empty and no terms appear on 
the right sides of those equations.) 

Here, each equation is solved for a principal derivative of maximal class k in such a way 
that the corresponding right side of the equation may depend on an arbitrary subset of the 
independent variables, an arbitrary subset of the dependent variables with 1 < /3 < Po, those 
derivatives uj for all 1 < 7 < m which are of a class j < k and those derivatives which are of the 
same class k but are not principal derivatives. Note that a principle derivative u'^ may depend 
on another principle derivative uj as long as ^ < k. The equations are grouped according to 
their class in descending order. 

Theorem 1 (Cartan-Kahler). Let the involutive differential equation TZi be locally represented 
in 5-regular coordinates by the system (f8a|) . ([Sbll . (|8c|l . Assume that the following initial condi- 
tions are given: 

n"(x\...,2:") = r(x\...,x"), /3i"^<a<m; (9a) 

n"(x\...,x"-\0) = r(x\...,x"-i), Z?}""^^ <a</3{"^; (9b) 
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u"{x\0,...,0)=r{x') 
n"(0,...,0)=r, 



/sf^ < a < /?; 

1 < a < 



(9c) 
(9d) 



'1 ' 



// the functions (j)'^ and are (real-) analytic at the origin, then this system has one and only 
one solution that is analytic at the origin and satisfies the initial conditions 

Proof. For the proof, see [19] or |23l Section 9.4] and references therein. The strategy is to 
split the system into subsystems according to the classes of the equations in it (see below) . The 
solution is constructed step by step; each step renders a normal system in fewer independent 
variables to which the Cauchy-Kovalevskaya theorem is applied. Finally, the condition that TZi 
is involutive leads to further normal systems ensuring that the constructed functions are indeed 
solutions of the full system with respect to all independent variables. ■ 

Under some mild regularity assumptions the algebraic equations can always be solved locally. 
From now on, we will assume that any present algebraic equation has been explicitly solved, 
reducing thus the number of dependent variables. We simplify the Cartan normal form of 
a differential equation as given in Definition [7] into the reduced Cartan normal form. It arises by 
solving each equation for a derivative u", the principal derivative, and eliminating this derivative 
from all other equations. Again, the principal derivatives are chosen in such a manner that their 
classes are as great as possible. Now no principal derivative appears on a right side of an equation 
(whereas this was possible with the non-reduced Cartan normal form of Definition [Tj) . All the 
remaining, non-principal, derivatives are called parametric. Ordering the obtained equations by 
their class, we again can decompose them into subsystems: 



Note that the values /?} are exactly those appearing in the Cartan test d?]), as the symbol 
matrix of a differential equation in Cartan normal form is automatically triangular with the 
principal derivatives as pivots. 

(k) (k) 

Definition 8. The Cartan characters of TZi are defined as a\ = m — PI and thus equal the 
number of parametric derivatives of class k and order 1. 

Provided that (5-regular coordinates are chosen, it is possible to perform a closed form invo- 
lution analysis for a differential equation TZi in reduced Cartan normal form. We remark that 
an effective test of involution proceeds as follows (see for example \23\ Remark 7.2.10]). Each 
equation in (jlOp is prolonged with respect to each of its non-multiplicative variables. The arising 
second-order equations are simplified modulo the original system and the prolongations with 
respect to the multiplicative variables. The symbol A/i is involutive if, and only if, after the 
simplification none of the prolonged equations is of second-order any more. The equation TZi is 
involutive if, and only if, all new equations even simplify to zero, as any remaining first-order 
equation would be an integrability condition. 

In order to apply this test, we now prove two helpful lemmata. We introduce the set B := 
uf is a principal derivative), and for each {a, i) B we define := — (p^- 
Using the contact fields ([T]), any prolongation of some (pf can be expressed in the following form. 

Lemma 1. Let the differential equation TZi be represented in the reduced Cartan normal form 
given by equation (|lUp . Then for any {a, i) B and 1 < j < n, we have 




< 



I < j < k < n, 
1 < a < /3f \ 
Pi^ < 7 < m. 



(10) 



m 
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Proof. By straightforward calculation; see [Gj Lemma 2.5.5]. 



For j > i, the prolongation Dj^f is non-multiplicative, otherwise it is multiplicative. Now 
let j > i, so that equation (jlip shows a non- multiplicative prolongation, and assume that we are 
using ^-regular coordinates. According to our test, the symbol A/i is involutive if, and only if, 
it is possible to eliminate on the right hand side of (jlip all second-order derivatives by adding 
multiplicative prolongations. 

If the differential equation is not involutive, then the difference 



(i) 



does not necessarily vanish but may yield an obstruction to involution for any (a, i) (z B and 
any i < j < n. The next lemma gives all these obstructions to involution for a first-order system 
in reduced Cartan normal form. 

Lemma 2. Assume that 5-regular coordinates are used. For an equation in Cartan normal form 
and indices i < j and a such that (a, i) G B, we have the equality 



a 



c 



E 



0\ 



i—1 m 

E E 



11^ 



J2 c!;{4>t)c^{<p]) 



(12) 



(13) 



PI 



(fe) 



E E 

l<h<k<i I 



^hk 



J2 c^,mc^{<t^i 



+ E ""hk 



(14a) 



(14b) 



i-l 

E 

h=l 



(0 



E 



U) 



+ E 



4^ 



/3P I3\-' 



7=/?r'+l 



7=/3l'''+l 



(15a) 



(15b) 



i-l 



E E 

i+l<k<j ''-f^l + 



^hk 



U) 



E 



7=/3f^+l 



(16) 
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j-l m 

-E E 

k=i s=fs[>'^+i 
i— 1 m 

-E E 



E 



0) 



(i) 



(17) 



(18) 



(19) 



Proof. By a tedious but fairly straightforward computation, see [Gj pages 48-55]. 



In line ()12p we have collected all terms which are of lower than second order. Furthermore, 
none of the appearing second-order derivatives is of a form that it could be eliminated by 
adding some multiplicative prolongation. Hence, under our assumption of (5-regular coordinates, 
the symbol TVi is involutive if, and only if, all the expressions in square brackets vanish. The 
differential equation TZi is involutive if, and only if, in addition line (jl2p vanishes, as it represents 
an integrability condition. Thus Lemma [2] gives us all obstructions to involution for TZi in 
explicit form. They will reappear in the proof of the existence theorem for integral distributions 
in Section [3 



5 The Vessiot distribution 



By Proposition [H the tangent spaces Tp{imjqa) of prolonged sections at points p G Jqir 
are always subspaces of the contact distribution Cq\p. If the section cr is a solution of the 
differential equation TZq, then by definition it furthermore satisfies im jgcr C TZq, and therefore 
T(im jgfj) C TTZq. Hence, the following construction suggests itself. 

Definition 9. The Vessiot distribution of a differential equation TZq C Jqir is the distribution 
V[7^g] C TTZq defined by 

Ti{V[nq]) =Ti{Tnq) nCq\n,. 

Note that the Vessiot distribution is not necessarily of constant rank along TZq (just like 
the symbol Mq); for simplicity, we will almost always assume here that this is the case. This 
definition of the Vessiot distribution is not the one usually found in the literature. But the 
equivalence to the standard approach is an elementary exercise in computing with pull-backs. 

Proposition 3. The Vessiot distribution satisfies V[TZq] = {i*Cq)^. 



JqTT is 



For a differential equation given in explicitly solved form, the inclusion map l: TZq 
available in closed form and can be used to calculate the pull-back of the contact forms. This 
has the advantage of keeping the calculations within a space of smaller dimension, namely the 
submanifold TZq. Thereby regarding the differential equation as a manifold in its own right, we 
bar its coordinates to distinguish them from those of the jet bundle. 

Example 1 . Consider the first-order system given by the representation TZi : Uf — v = Vf — Wx = 
Ux — w = 0. Then from the prolongations Uxt — Vx = and Uxt — wt = follows the integrability 
condition wt = Vx by elimination of the second-order derivative Uxt- Thus, the first projection 
of TZi is (in Cartan normal form) represented by 



TZ?: 



ut 

Ux 



- w. 



Vt 



Wx 



Wt 
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It is not difficult to verify that the projected equation 7^^^"* is involutive. For coordinates on Jitt 
choose X, t; u, v, w; Ux, Vx, Wx, ut, vt, wt- Since IZ^^^ is represented by a system in solved form, 
it is natural to choose appropriate local coordinates for 7^^^\ which we bar to distinguish them: 
X, t] u, V, W; v^, W^. The contact codistribution for Jivr is generated by 

uj^ = du — Uxdx — utdt, oJ^ = dv — Vxdx — vtdt, uj^ = dw — Wxdx — wtdt. 

The tangent space TIZ^^^ is spanned by dx-, du, d^, c?u7, d^^ and Tl(T71^^^) therefore by 
the fields dx, dt, du, dy + dut , + du^ , dy^ + and dy,^ + dy^ . This space is annihilated by 

= dUx — dw, oj^ = dvx — dwt, uj^ = dwx — dvt, lo^ = dut — dv. 

These seven one-forms lo^, . . . ,lJ^ annihilate the Vessiot distribution V[7^^^^], which is spanned 
by the four vector fields 

Xx=dx^ Uxdu + Vxdy + Wxdy, + Vxdy, + Wxdu^ , = dy^ + dy,^ , 

X2 = dt + utdy + vtdy + wtdy, + vtdut + wtdy^ , X4 = dy^ + dyj^ . 

In local coordinates on TZ^i \ these four vector fields become 

X2 = dt + vdu + w^dy + v^dw, Xa, = dyj^. 

They satisfy = Xj, as a simple calculation using the Jacobian matrix for Tl shows. The 
vector fields Xi are annihilated by the pull-backs of the contact forms, l*uj^ = du — u^dx — uldi, 
dtp' = dv — Ix^dx — xtdi and = dw — w^dx — wldi (the pullbacks of the remaining four 

one-forms uo^, . . . ,uj^ trivially vanish on 7^["^^). 

For a fully nonlinear differential equation TZq, in particular an implicit one, this approach 
to compute its Vessiot distribution V[7^g] via a pull-back is in general not effectively feasible. 
However, applying directly our definition of V[7^q], it is easily possible even for such equations to 
determine effectively Ti{y\JZq]), in other words: to realize it as a subbundle of T{JqTi)\'ji^. The 
contact fields ([T]) form a basis for Cq. It follows that for any vector field X G V[7^q], coefficients 
a*, 6^ S J-{JZq), where 1 < i < n, 1 < a < m and = q, exist such that 

i,X = a'cf^ + (20) 

If the differential equation TZq is locally represented by = 0, where 1 < r < i, it follows from 
the tangency of the vector fields in V[7?.g] that d^'^i^i^X^ = l^:X{^'^) = and thus the coefficient 
functions must satisfy the following system of linear equations: 

c('')(<J-)a* + C^(<J-)6^ = 0, (21) 

where 1 < r < i. Note that this approach to determine the Vessiot distribution is closely re- 
lated to prolonging the differential equation TZq and requires essentially the same computations. 
Indeed, the formal derivative ([5]) can be written in the form 

= C('')(<I>-) + C^(c^O^^+i, = (22) 

and in the context of the order-by-order determination of formal power series solutions (see 
for example |23^ Section 2.3]) these equations are considered an inhomogeneous system for the 
Taylor coefficients of order q + 1 depending on the lower order coefficients. Taking this point 
of view, we may call ()2ip the "projective" version of ()22p . In fact for n = 1, that is, for ordinary 
differential equations, this is even true in a rigorous sense. 
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Example 2. We consider the fully nonlinear first-order ordinary differential equation TZi locally 
defined by {u')"^ + v? + = 1. The contact distribution Ci is spanned by the two vector fields 
^1 = dx + u'du and X2 = du> and the Vessiot distribution Ti{y\JZq\) consists of all linear 
combinations of these two fields which are tangent to TZi. Setting lo = xdx + udu + u'du' , we 
thus have to solve the linear equation uj{aXi + 6X2) = in order to determine Ti{y\JZq\). Its 
solution requires a case distinction (which is typical for fully nonlinear differential equations). 
If u' 7^ 0, then we find 

Ti{V[nq]) = {u'dx + {u'fdu - (X + U'u)du'). 

For u' = and x ^ 0, the Vessiot distribution is spanned by the vertical contact field X2. 
Finally, for x = u' = the rank of the Vessiot distribution jumps, as at these points the whole 
contact plane is contained in it. 

The definition of the symbol A/'g and of the Vessiot distribution V[TZq], respectively, of a diffe- 
rential equation TZq C J^tt immediately imply the following generalisation of the above discussed 
splitting of the contact distribution Cq = Vtt^_-^ @'H (such a splitting of the Vessiot distribution 
is also discussed by Lychagin and Kruglikov |141 [T5] where the Vessiot distribution is called 
"Cartan distribution"). 

Proposition 4. For any differential equation TZq, its symbol is contained in the Vessiot distri- 
bution: Afq C V[7^g]. The Vessiot distribution can therefore be decomposed into a direct sum 

VlTZq] =Mq®n (23) 

for some complement Ti {which is not unique). 

Such a complement Tl is necessarily transversal to the fibration TZq — > X and thus leads 
naturally to connections: provided dim 7^ = n, it may be considered the horizontal bundle of 
a connection of the fibred manifold TZq — > X. 

Definition 10. Any such connection is called a Vessiot connection for TZq. 

In general, the Vessiot distribution V[7^q] is not involutive (that is, closed under the Lie 
bracket; an exception are formally integrable equations of finite type [23\ Remark 9.5.8]), but it 
may contain involutive subdistributions. If these are furthermore transversal (to the fibration 
TZq — > X) and of dimension n, then they define a flat Vessiot connection. 

Lemma 3. // the section a : X ^ £ is a solution of the equation TZq, then the tangent bundle 
T{imjqa) is an n-dimensional transversal involutive subdistribution ofV\TZq\\imjq(T- Conversely, 
ifUCl V[TZq] is an n-dimensional transversal involutive subdistribution, then any integral man- 
ifold ofU has locally the form imjgfj for a solution a of TZq. 

Proof. Let u be a local solution of the equation TZq. Then it satisfies by Definition [3] imj^u C 
TZq and thus T{imjqa) C TTZq. Besides, by the definition of the contact distribution, for x ^ X 
with jqcr(x) = p S Jgvr, the tangent space Tp{[mjqa) is a subspace of Cq\p. By definition of the 
Vessiot distribution, it follows Tp{imjqcr) C Ti{TpTZq) fl Cg|p, which proves the first claim. 

Now let U C V[7?.q] be an ?i-dimensional transversal involutive subdistribution. Then ac- 
cording to the Frobenius theorem, U has n-dimensional integral manifolds. By definition, 
Ti{y\JZq\) C Cq\Ti^; this characterises prolonged sections. Hence, for any integral manifold of ^ 
there is a local section a such that the integral manifold is of the form imjgcr. Furthermore, the 
integral manifold is a subset of TZq. Thus it corresponds to a local solution of TZq. ■ 

This simple observation forms the basis of Vessiot 's approach to the analysis of TZq: he 
proposed to construct all flat Vessiot connections. Before we do this, we first show how integral 
elements appear in this program. 
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Proposition 5. Let lA C V[7^g] he a transversal subdistribution of the Vessiot distribution of 
constant rank k. Then the spaces Up are k- dimensional integral elements for all points p € IZq 
if, and only if, [U,l/{] C V[7?-g]. 

Proof. Let {cJi, . . . ,uJr} be a basis of the codistribution L*Cq. Then an algebraic basis of the 
ideal I[TZq] is {u>i, . . . , ujr,du>i, . . . , duJr}- Any vector field X ^ U trivially satisfies u>i{X) = 
by Proposition [3l For arbitrary fields Xi,X2 € U, we have 

du;iiXi,X2) = XiicoiiX^)) - X2{uj^{Xi)) + LUi{[Xi, X2]) . 

The first two summands on the right hand side vanish trivially and the remaining equation 
implies our claim. ■ 

We call a subdistribution U C V[7^g] satisfying the conditions of Proposition O an integral 
distribution for the differential equation TZg. In the literature [2l], the terminology "involution" 
is common for such distributions which, however, is confusing. Note that generally an integral 
distribution is not integrable; the name only reflects the fact that it consists of integral elements. 

A general differential equations TZg does not necessarily possess any Vessiot connection (not 
even a non-flat one). Their existence is linked to the absence of integrability conditions. More 
precisely, we obtain the following characterisation. 

Proposition 6. Let TZg be a differential equation. Then its Vessiot distribution possesses locally 
a direct decomposition with an n-dimensional complement 7i such that V\TZq\ = Mq (BTi if, and 
only if, there are no integrability conditions which arise as the prolongation of equations of lower 
order in the system. 



Proof. Let TZg be locally represented by 



<^^(x,u(9)), 
ip"^fx,u(«-^)l 



such that the equations ^"^(x, u^'^^) = do not imply lower-order equations which are indepen- 
dent of the equations W^{x,u^''~^^) = 0. Let Uj-^) denote the subset of all derivatives of order q 
only; then the Jacobi matrix (9<?'^(x, u'^'?^)/(9u(g)) has maximal rank. If we proceed as in the 
last proof, then the ansatz (j20p for the determination of the Vessiot distribution yields for the 
above representation the linear system 

cl''\<P^)a' + C^{<P^)b'^ = 0, (24) 

Here, the matrix Ca{0'^) has maximal rank, too; thus the equations c'f\'^'^)a^ + Ca{^^)b'^ = 
can be solved for a subset of the unknowns 6^. And since no terms of order q are present in 

!P'°'(x, u'-'?" ') = 0, we have Cf''^(l^'^) = A?^'"- Recall that we consider the Vessiot distribution, 
and thus the linear system (j24|) . only on TZg. It follows that the subsystem c\'^^ {^")a^ = 
becomes trivial if, and only if, no integrability conditions arise from the prolongation of lower 
order equations. And if, and only if, this is the case, then (j24p has for each 1 < j < n a solution 
where = 1 while all other a* are zero. The existence of such a solution is equivalent to the 
existence of an n-dimensional transversal complement 7i. ■ 

Prom the proof of this proposition follows that for the determination of the Vessiot distri- 
bution V[7^g], equations of order less than q in the local representation of TZg can be ignored 
if there are no integrability conditions which arise from equations of lower order. It is the 
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integrability conditions which arise as prolongations of lower order equations that hinder the 
construction of n-dimensional complements, while those which follow from the relations between 
cross derivatives do not influence this approach. 

Remark 1. If one does not care about the distinction between different kinds of integrability 
conditions and simply requires that TZq = TZ^q^ (meaning that no integrability conditions at all 
appear in the first prolongation of 7tq), then one can provide a more geometric proof for the 
existence of an n-dimensional complement (of course, in contrast to Proposition [U the converse 
is not true then). 

The assumption TZq = TZ^q^ implies that to every point p G TZq at least one point p G Tlq+i 
with TTq~^^{p) = p exists. We choose such a point and consider imTq^ilp) C Tp{Jq7r). By 
definition of the contact map Tg+i, this is an n-dimensional transversal subset of Cq\p. Thus 
there only remains to show that it is also tangential to TZq, as then we can define a complement 
by Tpi{Hp) = \mTqj^i{p). But this tangency is a trivial consequence of p G using for 

example the local coordinates expression ([3]) for the contact map and a local representation 
= of TZq, one immediately sees that the vector Vi = Tqj^i{p,dxt) G Tp{JqTi) satisfies 
d(t'^\p{vi) = Di^'^{p) = and thus is tangential to TZq. 

Hence it is possible to construct for each point p G TZq a transversal complement Tip such that 
= iMq)p(BTi.p. There remains to show that these complements can be chosen so that they 
form a smooth distribution. Our assumption TZq = TZ^q^ implies that the restricted projection 
TTg"^^ : TZq^i TZq IS a surjective submersion, that is, it defines a fibred manifold. Thus if we 
choose a section 7 : TZq — > TZq+i and then always take p = j{p), it follows immediately that the 
corresponding complements Tip define a smooth distribution as required. 

Example 3. Consider again the differential equation TZi in Example [TJ It is locally represented 
by the same equations as TZ^i \ except that the integrability condition wt = Vx is missing. The 
matrix of Tl for the system TZi has eleven rows and eight columns — one column more than 
the symbol matrix for the system TZ^^\ The symbol Ti(7Vi) of TZi is 3-dimensional, spanned 

by dy^, dwt and dw^ + d^^., while the symbol Ti{J\fi^^) of TZ^^^ has dimension 2 and is spanned 
by dy.^ + and + d^^. But the one-forms w^, and lu'^ (and their pull-backs i*uj^ = 
du — u^dx — ut(R, CuP' = dv — v^dx — vtdt and l*uj^ = dw — w^dx — wldi) are the same, and 
therefore the coordinate expressions for the vector fields Xi and X2 (and their representations 

Xi = dx + u^du + v^dv + vJ^du) and ^2 = 9^ + utdu + vldv + wtdw in TTZi and TTZ^^ ) look alike 
(see Example [1] for their representations). The integrability condition wt = Vx does not influence 
the results as it stems from the equality of the cross derivatives, utx = Vx and Uxt = wt, and not 
from the prolongation of a lower order equation. 

Now consider for comparison the differential equation which is locally represented by 

{Ut = V, Vt = Wx, Wt = Vx, 

Ux = W, 
U = X. 

It arises from the system 7^^^^ in Example [T] by adding the algebraic equation u = x. Proceeding 
as in Example [H we find that the Vessiot distribution V[7?-i] is spanned by the three vector 
fields 

d^, cfer and vdx + {I -w)di+ {uj^{l -w) + vv^^ 

The vector fields and d^^ generate the symbol Mi; any of its complements in V[7^i] is one- 
dimensional and, as the dimension of the base space is two, none of them has the right dimension 
to be the horizontal space of a connection. 
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The reason for this is that TZi is not formahy integrable, as the prolongation of the algebraic 
equation u = x leads to the integrability conditions Ux = 1 and ut = 0. Projecting the prolonged 
equation gives 

{Ut = V = Vt = Wx = Wt = Vx = 0, 
Ux = w = 1, 
U = X. 

Now the symbol vanishes, and so do the pull-backs of the contact forms: i*uj^ = dx — dx = 0, 
L*u;^ = dv = 0, i*a;3 = dw = 0. Therefore we find VlTZ^^^] = Af^^^ e 7^ = {0} {ds, d^). As the 

Lie brackets of dx and trivially vanish, TTZ^^^ = V[iZ^i^] =1-1= {dx,di) is a two-dimensional 
involutive distribution. 

Any ?i-dimensional complement Ti is obviously a transversal subdistribution of V[7?-g], but not 
necessarily involutive. Conversely, any n-dimensional subdistribution Ti of V[7^g] is a possible 
choice as a complement. Choosing a "reference" complement TIq with a basis (Xj : 1 <i<n), 
a basis for any other complement 7i arises by adding some vertical fields to the vectors Xi. 
We will follow this approach in the next section. For the remainder of this section we turn our 
attention to the choice of a convenient basis of V[7^g] that will facilitate our computations. 

Let r := dimTVg. As an intersection of two involutive distributions, the symbol Afq is an 
involutive distribution, too. Hence, there exists a basis (Y^ : 1 < A: < r) for it whose Lie brackets 
vanish: [^^,1^] = for all 1 < k,£ < r. Since the vertical bundle V7r^_^ is also involutive, we 
can decompose it into 

where W is again an involutive distribution. It can be spanned by vector fields Wi, W2, ■ ■ ■ , Ws 
(where s = "^^=1 M'^^ equals the number of principal derivatives) which are chosen such that 
we have [Wa, W^] = for all 1 < a, 6 < s. In local coordinates, a particularly convenient choice 
for the fields Yfc and Wa exists. We first choose for any 1 < k < r a parametric derivative 
that is (a,/i) ^ B, and set Yfc := := t,,(c^); then we choose for any 1 < a < s a principal 
derivative that is {a,n) € B, and set Wa ■= W^ := du'^. 

The reference complement TIq is chosen as follows. Any basis of it must consist of n transversal 
contact fields. Since the fields Co are vertical, we can always use a basis {Xi : \ <i <n) oi the 
form 

x, = c\'^^ + if^c^ 

with some coefficient functions ^f^ chosen such that Xi is tangential to IZq. The fields Ca also 
span the vertical bundle Vtt^_i and hence we may exploit the above decomposition for a further 
simplification of the basis. By subtracting from each Xi a suitable linear combination of the 
fields Yfe spanning the symbol Mq, we arrive at a basis (X,: 1 <i <n) where 

X, = Ci"^ + CtWa. (25) 

As already mentioned above, the Vessiot distribution V[7?.g] is not necessarily involutive. So 
it is not surprising that its structure equations are going to be important later. We may extend 
the above chosen basis (Xj, Y^) of V[7^g] to a basis of the derived Vessiot distribution, V'[TZq], 
by adding vector fields Zc, 1 < c < C := dimV'[7^g] — dimV[7^g], where, using ([2|), for each c 
we have some coefficients such that Zc = K^pdu'^ with \v\ = q — 1. By construction, the 
non- vanishing structure equations of V[7^q] take now the form 



[Xi,Xj] = ef^Zc and [X^,Yk] = Hf^Z, 



(26) 
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for 1 < i,j < n and 1 < /c < r, with smooth functions O^- and S'^^. (For the complete set of 
structure equations, we have to add [Y^, YJ] = for 1 < k^l < r.) 

Remark 2. Since the vector fields Z^., which appear on the right sides of the structure equa- 
tions (j26p . may, for q = 1, span a proper subspace in (9„q : 1 < a < m), about the exact size 
of which we know nothing, we write them as linear combinations Z^. =: K'^du^. The structure 
equations then become 

[Xi,Xj] = =: 0° , 1 < i, i < n, m^') 

[Xi, Yk] = =: ^f^a,. , 1 < i < n, l<k<r. m^') 

Knowing the larger sets of coefficients 0^-, 5"^^, we can reconstruct the true structure 
coefficients Ofj, S^f^ by solving the overdetermined systems of linear equations 

This is always possible since the fields Zc are assumed to be part of a basis for the derived 
Vessiot distribution V'[7^i] and therefore linearly independent. Thus there exist some coefficient 
functions such that 

For our later proof of an existence theorem for integral distributions, we will have to analyse 
certain matrices with the coefficients Ofj and S"?^ as their entries. It turns out that this analysis 
becomes simpler, if we use the extended sets of coefficients Ofj and S^f^ instead. 

For a first-order equation TZi with Cartan normal form ^ satisfying the assumptions of 
Proposition [6] it is possible to perform this process explicitly. We choose as a reference comple- 
ment Tlo the linear span of the vector fields 

..X. = C(^)+ C^'\<P'',)Cli (27) 

for 1 < i < n. One verifies in a straightforward computation that (j27p represents a valid choice 
(see [H Proposition 3.1.19]). Using this reference complement, we can explicitly evaluate the Lie 
brackets (j26p on TZi . As we are not able to determine a simple expression for the derived Vessiot 
distribution, we follow the approach of Remark [2] and consider the equations (j26l ) instead. 

Lemma 4. Let i < j, without loss of generality. Then we obtain for the extended set of structure 
coefficients 0^- in local coordinates on TZi the following results: 

Ci'\<pj) - Cf\(t>t), («,^) G B and (a, j) e B, 

C[^\<P';), (a,i) B and {a,j) G B, (28) 

0, {a,i) ^ B and {a,j) ^ B. 

Proof. By a rather straightforward computation; see [U pages 75 and 76]. ■ 

We collect these coefficients into vectors Oij which have m rows each where the entries are 
ordered according to increasing a. It is useful to denote the symbol fields Yfe = i^,{d-j) by using 

double indices: Ife = vj^ for any (/?, h) B. 



0" = < 
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Lemma 5. Set =: Yf^ , then the extended set of structure coefficients S^j^ in local coordinates 
on TZi are 



0, 



(a, i) G B, 

{oi,i) ^ B and (a,i) = {(3,h), 
{a,i) ^ B and (a,z) 7^ {(3,h). 



(29) 



Proof. Again by a rather straightforward computation; see [6l page 76]. 



Example 4. We calculate the structure equations for Example [TJ Here, = 1 and = 3. 
We set X3 =:Yi = Y^ and X4 =: I2 = ^i'^- Then besides the trivial structure equations we get: 



?(2) 



[XuX2] = 0, [Xi,Yi] = [X2,Y2] 



-dv, 



[X2,Y^] = [X^,Y2] = -d^. 



Noting that here the set B contains the pairs (1, 1) = Ux, (1, 2) = ut, (2, 2) = vt and (3, 2) = wt 
while (2, 1) = Vx and (3, 1) = Wx are not in B, one easily verifies that all coefficients are as 
given in equations ([25]) and 

We end this section with two technical remarks on how these coefficients are collected into 
matrices S^. The examination of the ranks of these matrices Si is basic for the proof of the 
existence Theorem [2] for integral distributions. 



Remark 3. Some of the terms —C^{<j)f) where (a,i) S B vanish, too: all of the parametric 
derivatives on the right side of an equation <Pf in the reduced Cartan normal form (fTOj) are of 
a class lower than that of the equation's left side as otherwise we would solve this equation for 



the derivative of highest class. This means 



whenever i = cls{uf) < cls(uf ) = h. 



We collect the coefficients S^j^ into matrices Si using i as the number of the matrix to which 
the entry belongs, a as the row index of the entry and k as its column index. These matrices 
have m rows each, ordered according to increasing a; and they have r = dimTVi columns each 



of which can be labelled by pairs (/?, h) ^ B or the symbol fields Y^ 
1 < h < n, we set 



yP . More precisely, for 



(h) 



-c 



h 



(h) 



(0 



-a 



V 



(30) 



Such a matrix with an upper index h collects all those into a block where (a, i) € B. For 



any 1 < i < n, we have m — I3\ 



\ so such a matrix has rows and q^*^ columns. Since, 



for any h with i < h and for all P^'' + 1 < /? < m, we have —C^{(pf) 



?{^) , 1 ^ « ^ _ r<h^Aa\ _ g^^^.^^ matrices [Si]^ 

where i < h are zero. Furthermore, for 1 < h < n, we assemble the remaining terms S^j^ (which 
are those where {a,i) ^ B) in a matrix. As above, let Y^ = Yj^ , and denote any entry S°j^ 
accordingly, for the sake of introducing the following matrix, by ^S". Now set 



^/^r'+i h 



(i) 



h ^l^i 



h 



h 

m'- 



7^1 



(i) 



h 
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(0 



„(') 



~a 



(2) 



Figure 1. A sketch for the matrix in equation PT|) . 



For any 1 < h < n, such a matrix with the index h written below has a^*^ rows and a^^'^ 
columns. Let for any natural numbers a and h denote 0^x6 the o x 6 zero matrix. According to 
equation (j29p . we have 



-1 



[■='i\h 



(0 V, Ci) ) 



h = i, 



Using the matrices [Si]^ and [Si]h as blocks, we now build the matrix 



2 



Taking into account what we have noted on its entries, this means 








-1 



0---0 



(31) 



A sketch of such a matrix Si where the entries which may be different from zero are marked as 
shaded areas and — 1 (i) as a diagonal line is given in Fig. [H 



For all h where 1 < h < n, we call the block [S'j] in Si stacked upon the block [Si]h in Si 
the hth block of columns in Si. For those h with /jj^^ = m the hth. block of columns is empty. 
Now the symbol fields Y"^, or equivalently the pairs {f3, h) ^ B, are used to order the dimA/i = r 



columns of Si, according to increasing h into n blocks (empty for those h with a\ 
within each block according to increasing 13 (with + 1 < ,9 < m). 



(h) 



0) and 



Remark 4. This means, the columns in Si are ordered increasingly with respect to the term- 
over-position lift of the degree reverse lexicographic ranking. Therefore, the entry — C'^^^j {4>i) 

stands in the matrix Si in line a, in the hth block of columns of which it is the 7th one from 
the left. Entries different from zero and from —1 may appear in Si only in a [Si]^ for h < i. 

" (i) 

To be exact, for any class i, the matrix Si has a\ rows where all entries are zero with only one 
exception: for each 1 < ii < a^^ we have 
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where I := '^h=i'^i^^ + "^^^ entries in the remaining upper rows are —C^{(j)'f). The 
potentiahy non-trivial ones of them are marked as shaded areas in Fig. [TJ 

Note that for a differential equation with constant coefficients all vectors Oij vanish and for 
a maximally over-determined equation there are no matrices S^. 

(i) 

The unit block of a\ rows, 1^(0; leads immediately to the estimate 




< rank Si < min < m, 



Example 5. For Example [U we have n = 2 and therefore two matrices Si and S2. Their 
entries follow immediately from our results in Example SI For i = 1, we have 



"1 




The first line is [Si]^ = (0, 0). The unit block below it is [S*!]!- (We have f3^^ = 1 and therefore 
a|^^^ = 2. There is only /i = 1 to consider, as both parametric derivatives Vx and Wx are with 
respect to x only.) Because of the very simple nature of our system, we find here by accident that 
Si = S2. However, [^2]^ = "2 is the whole matrix while there is no [^2)1 because (3^^ = m = 3 
and therefore a^i^ = 0. Finally, we obtain ©12 = {vx — Wf,0, 0)*. (The t top right marks the 
transpose.) Note that its first entry is in fact the integrability condition. 



6 Flat Vessiot connections 



In this section, we develop an approach for constructing flat Vessiot connections which improves 
recent approaches [SllMlllZ]: we exploit the splitting of V[7?.g] suggested by (p3]) to introduce 
convenient bases for integral distributions which yield structure equations that are particularly 
simple. Finally, we give necessary and sufficient conditions for Vessiot's approach to succeed. 

Let TZq locally be represented by the system 'P'^ {x,u^''^) = where 1 < r < t. Our goal 
is the construction of all n-dimensional transversal involutive subdistributions U within the 
Vessiot distribution V[7^g]. Taking for the Vessiot distribution the basis (Xj,Yfc) where the 
vector fields are the above mentioned basis of the symbol TL[Mq) with vanishing Lie brackets 
and the vector fields Xi are given in (|25p . we make for the basis ([/, : 1 < i < n) of such 
a sub distribution U the ansatz 

Ui = X, + C^Yk 

with yet undetermined coefficient functions G ^{T^q)- This ansatz follows naturally from 
our considerations so far, as the fields Xi are transversal to the fibration over X and span a 
reference complement to the symbol and in Afq all fields Y^ are vertical. 

Since the fields Ui are in triangular form, the distribution U is involutive if, and only if, their 
Lie brackets vanish, and using the structure equations (j26|) this means: 

[u,,Uj] = [x,,x,] + c^[Yk,Xj] + c^[x„Y,] + {u,ic^) - Ujich)yk 

= (0f,. - S^,C^ + S^C^)Z, + (f/,(Cj=) - U,{Ch)Y, = 0. (32) 

It follows from the definition of the fields and Zc that they are linearly independent and so 
their coefficients must vanish for U to be involutive. Thus the Lie bracket ([32]) yields two sets 
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of conditions for the coefficient functions a system of algebraic equations 



and a system of differential equations 



H!^:=U.iq)-U,iCn = 0, { (34) 

In the algebraic system P5|) the true structure coefficients Ofj, appear. For our subsequent 
analysis we follow Remark [2] and replace them by the extended set of coefficients ©j", SJ^. This 
corresponds to replacing (j33|) by an equivalent but larger linear system of equations which is, 
however, simpler to analyse. 

Obviously, (j33|) is an inhomogeneous linear system to which any solution method for linear 
systems can be applied. We shall see in the next section that its structure (induced by the 
structure equations of the Vessiot distribution) allows us to decompose it into simpler subsystems 
which are considered step by step. Many papers on the Vessiot theory (see for example [5} 124]) 
study at this stage a quadratic system which only by such a step-by-step approach can be 
reduced to a series of linear problems. The linearity of (jSSp is a simple consequence of our choice 
of a basis for V[7^g] which in turn exploits the natural splitting V[7?-q] = Mq © H. 

Remark 5. The vector fields lie in the Vessiot distribution V[7^i]. Thus, according to 
Proposition [51 the subdistribution U is an integral distribution if, and only if, the coefficients Q 
satisfy the algebraic conditions (jSSp . This observation permits us immediately to reduce the 
number of unknowns in our ansatz. Assume that we have values 1 < i,j < n and 1 < a < m 
such that both {a,i) and (a,j) are not contained in B, that is, uf and uj are both parametric 
derivatives (and thus obviously the second-order derivative ufj, too). Then there exist two 
symbol fields = i*{du'^) and Yi = i*(c^u'?)- Now it follows from the coordinate form ([3|) of 
the contact map that the subdistribution U, spanned by the vector fields Uh = Xh + Ch^k for 
1 < h < n, can be an integral distribution if, and only if, 

Cj= = d or, equivalents, c]°"^ = Ci°''^ (35) 

for all 1 < i < j < n and 1 < k,l < r. 

As the unknowns Cj ™ay be understood as labels for the columns of the matrices S^, this 
identification leads to a contraction of these matrices. The contracted matrices, denoted by 
Sh, arise as follows: whenever ^j^ = then the corresponding columns of are added; their 
sum replaces the first of these columns, while the second column is dropped. Similarly, we 
introduce reduced vectors where the redundant components are left out. From now on, we 
always understand that in the equations above this reduction has been performed. Otherwise 
some results would not be correct (see Example 1101 below) . 

For a more thorough outline of these technical details, see [6j Lemma 3.3.5 and Remark 3.3.6]. 

Example 6. We consider a first-order equation in one dependent variable 

Un = </'n(x,n, Ui, . . .,Un-r-l), 

7^l: ; 

'^n—r — 4'n—r (^j "U, 1*1 , . . . , l*n— r— l) ■ 



It is easy to show that the symbol of such an equation is always involutive and that the used 
coordinates are 5-regular. In order to simplify the notation, we use the following convention for 
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the indices: 1 < k,i < n — r and n — r < a,b < n. Thus we may express our system in the 
concise form Ua = <j)a{^,u,Uk) and local coordinates on the submanifold TZi are (x, U, H^). 

The pull-back of the contact form uj = du — Uidx'^ generating the contact codistribution is 
= du — Ukdx^ — (padx"" and V[7^i] = {Xi, . . . , Xn,Yi, ■ ■ ■ , ^n-r-i) where 

Xk = d^k + Ukdu, Xa = dx'^ + <Padu, % = ■ 

The fields span the symbol A/i and the fields Xi our choice of a reference complement Ti.Q. 
Setting Z = du, the structure equations of V[7^i] are 

[Xk,Xi]=0, [Yk,Ye]=0, 

[Xk,Xa] = Xki^a)Z, [Xa,Xh] = {XaiM " M^a)) Z , 

[Xk,Ye] = -5mZ, [Xa, Yk] = -Yk{(t>a)Z. 

Now we make the above discussed ansatz C/j = Xj + Q^Yk for the generators of a transversal 
complement Ti. Modulo the Vessiot distribution V[7^i] we obtain for their Lie brackets 

[Uk,U,] = {Ci-C!l)Z modV[7^l], (36a) 

Pa, Uk] = (Ca - Qiyi{<l>a) " Xk{(l>a)) Z mod VilZi], (36b) 

[Ua, f/fe] = (Can(<^6) - QiYi{(^a) + X^ - X^) Z mod V[7^l]. (36c) 

The algebraic system (j33p is now obtained by requiring that all the expressions in parentheses 
on the right hand sides vanish. Its solution is straightforward. The first subsystem (I36ap implies 
the equalities Ci ~ Ce- This result was to be expected by the discussion in Remark [5} both Uk 
and ue are parametric derivatives for TZi and thus we could have made this identification already 
in our ansatz for the complement. The second subsystem ()36bp yields that = C|^(0a) + 
Xki4>a)- If we enter these results into the third subsystem ()36cp . then all unknowns C drop out 
and the solvability condition 

XaicPb) - XM + Xk{cPa)Yk{<Pb) - Xk{<Pb)Yk{(l>a) = 

arises. Thus in this example the algebraic system (j33p has a solution if, and only if, this condition 
is satisfied. 

Comparing with the classical theory of such systems, one easily verifies that this solvability 
condition is equivalent to the vanishing of the Mayer or Jacobi bracket [ua — (j)a,Ub — (pb] on the 
submanifold TZi, which in turn is a necessary and sufficient condition for the formal integrability 
of the differential equation TZi. Thus we may conclude that TZi possesses n-dimensional integral 
distributions if, and only if, it is formally integrable (which in our case is also equivalent to TZi 
being involutive). 

Thus here involution can be decided solely on the basis of the algebraic system (j33p . We will 
show in the next section that this observation does not represent a special property of a very 
particular class of differential equations but a general feature of the Vessiot theory. 



7 The existence theorem for integral distributions 

Now the question arises, when the combined system ()33p . (j34p has solutions? We begin by 
analysing the algebraic part (f33l) . We use for this analysis a step-by-step approach originally 
proposed by Vessiot [28]. Our analysis will automatically reveal necessary and sufficient assump- 
tions for it to succeed. As outlined in Remark [5l we replace C2 by C2 since for the entries ci^'^^ 
where + 1<P< m we know already from equation ([35|) that ci^'^^ = Cj^'^^- Thus we begin 
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the construction of the integral distribution lA by first choosing an arbitrary vector field Ui and 
then aiming for another vector field U2 such that [f/i, Uq] € Ti{y[JZq\). During the construction 
of the field U2 we regard the components of the vector Qi = C,i as given parameters and the 
components of C2 as the only unknowns of the system 



^lC2 = ^2Cl - O 



12- 



(37) 



Since the components of Qi are not considered unknowns, the system ([37]) must not lead to any 
restrictions for the coefficients Ql. Obviously, this is the case if, and only if, 

rank ^1 = rank (Si S2). (38) 

Assuming that ([38|) holds, the system ([37|) is solvable if, and only if, it satisfies the augmented 
rank condition 



rank Si = rank (Si S2 —Oi2)- 



(39) 



When we have succeeded in constructing U2-, the next step is to seek a further vector field U-j, 
such that [[/i,f/3] € Ti(V[7^g]) and [C/2,C/3] G Ti{y[R.q\). Now the components of both vectors 
Ci and C2 are regarded as given, and the components of (^3 are regarded as the unknowns of the 
system 



S1C3 = S3C1 - e 



13 1 



"2C3 — "3C2 - & 



23- 



Now this system is not to restrict the components of both C,i and C,2 any further. The interrela- 
tions between the Q following from the condition (j35p for the existence of integral distributions, 
given in Remark [U are taken care of by contracting S3 into S3. This implies that now the rank 
condition 



rank 



"1 

"2 



rank 



Si S3 
^2 S3 



has to be satisfied. If it is, then for 1 < c < C = dim V'[7?-g] — dim V[7^q] the system 



e 



13 



"3fcSl -r "lfc'53 



0, 



0: 



23 



a/c I '^c /-k 
"3fcS2 "T "2fcS3 



is solvable if, and only if, the augmented rank condition 



rank I ^ ] = rank 



Si S3 
^2 S3 



-013 
-023 



holds. Now we proceed by iteration. Given j — 1 vector fields Ui, U2, ■ ■ ■ , Uj-i of the required 
form spanning an involutive subdistribution of Ti(V[7^i]), we construct the next vector field Uj 
by solving the system 



SiO = sXi-0 



(40) 







Again we consider only the components of the vector Q unknowns, and the system (|40p must 
not imply any further restrictions on the components of the vectors Ci for 1 < i < j. The 
corresponding rank condition is 



rank 



/ ^1 \ 

"2 

V^.-i/ 



rank 



/ Ji 

"2 



\ 



s,7 



(41) 
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Assuming that it holds, the equations (j40p are solvable and yield solutions for the components 
of Cj if, and only if, the system satisfies the augmented rank condition 



rank 



rank 



: 



-02, 



(42) 



u 



Remark 6. We can simplify our calculations for a differential equation TZq, represented by 
> — '/'^(X' ^1 u^'^^) where each equation is solved for a principal derivative (all pairwise different) 
u'j^ with I /Li I = q, and where u^'^-* denotes the set of the remaining, thus parametric, derivatives 
of order q and less. (Example [U discusses a system of that kind.) Now we can use the local 
coordinates on TZg. For the generators of the symbol Mq we may choose for (a, n) ^ B the vector 
fields = as a basis for the complement H C V[7?.q], we can take the vector fields 



a=l 0<\fi\<q 



d,,a , 



which satisfy equation (I27p . For the vector fields Wa which appear in equation (j25p we may 
choose the contact vector fields Ca where (a, //) G B. The further procedure, solving first the 
structure equations, which now take the simple form ([26]) . to find the generators of V'iJZq) and 
then solving equation ()33p for the coefficient functions Cjj is the same as in the case where the 
representation is not in the above solved form. 

Example 7. Before we proceed with our theoretical analysis, let us demonstrate with a concrete 
differential equation that in general we cannot expect that the above outlined step-by-step 
construction of integral distributions works. In order to keep the size of the example reasonably 
small, we use the second-order equatioi|l] 7^2 defined by the system Uxx = c^u and Uyy = fiu 
with two real constants a, (3. Its symbol M2 is not involutive. However, one easily proves that 
7^2 is formally integrable for arbitrary choices of the constants a, [3. One readily computes that 
the Vessiot distribution V[7^2] is generated by the following three vector fields on TZ2'- 

Xi = dx + u^du + audu^ + Th^du^, X2 = dy + u^du + u^du^ + Yi = du^. 

They yield as structure equations for V[7?.2]: 

[Xi,X2] = IM^du^ - au^du^, [Xi,Yi] = [X2, Yi] = -d^^. 

For the construction of a two-dimensional integral distribution U C V[7^2] we make as above 
the ansatz Ui = Xi + C,}Yi with two coefficients Q. As we want to perform a step-by-step 
construction, we assume that we have chosen some fixed value for C,\ and try now to determine C,\ 
such that [C/i, U2\ = mod V[7?.2]. Evaluation of the Lie bracket yields the equation 



y 



(43) 

Hence one cannot 



[f/i, U2] = {fiui- C2)% - («% - C,l)du^ mod V[7^2]. 
A necessary condition for the vanishing of the right side is that C,\ 

choose this coefficient arbitrarily, as we assumed, but ()43p determines both functions C} uniquely. 
Note that the conditions on the coefficients C} imposed by ([13|) are trivially solvable and thus 7^2 
possesses an integral distribution (which is even involutive). The problem is that it could not 
be constructed systematically with the above outlined step-by-step process. 



'^It could always be rewritten as a first-order one satisfying the assumptions made above, and the phenomenon 
we want to discuss is independent of this transformation. 
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The following theorem links the satisfaction of the rank conditions (|4ip and (|42p . and thus 
the solvability of the algebraic system ([55]) by the above described step-by-step process, with 
intrinsic properties of the differential equation TZq and its symbol Mq. It represents an existence 
theorem for integral distributions. 

Theorem 2. Assume that 6-regular coordinates have been chosen for the local representation of 
the differential equation IZq. Then the rank condition (j4ip is satisfied for all 1 < j < n if, and 
only if, the symbol Mq is involutive. The augmented rank condition (l42]) holds for all 1 < j < n 
if, and only if, the differential equation TZq is involutive. 

The proof of Theorem [2] is given in Appendix |Al as parts of it are rather technical. First, we 
transform the differential equation into an equivalent first-order system with a representation in 
the reduced Cartan normal form (jlOp and the same Cartan characters. This is always possible; 
see, for example, [231 Proposition A.3.1]. The basic idea of the proof is then that the rank 
condition (j4ip is equivalent to the vanishing of the obstructions to involution of the symbol in 
Lemma O and the augmented rank condition (I42p is equivalent to the vanishing of the remaining 
integrability conditions (recall that Lemma [2] gives all obstructions to involution only if the used 
coordinates are 5-regular so that this assumption is necessary in Theorem [2]). The concrete 
implementation of this idea requires some technical considerations concerning the transformation 
of the matrices (j4ip and (j42p into row echelon form, working out their contractions and analysing 
the interrelation between these operations. 

Example 8. We demonstrate the role of ^-regularity of the coordinates with the wave equation 
Uxy = in characteristic coordinates which are not (5-regular. A straightforward computation 

(2) (2) 

yields as generators for the Vessiot distribution the fields Xi = C]; , X2 = C2 , Yi = du^^ 
and I2 = duyy- Following our approach, we make the ansatz Ui = Xi + ClYi + C1Y2 and 
U2 = X2 + C2Y1 + (2^2- Evaluation of the Lie bracket [Ui, U2] yields now the algebraic equations 
Ci^ = and C2 = and the differential equations Ui{C^) - C/2(Ci) = and Ui{Ci) - t^2(Ci) = 0- 
As in Example [71 we see that the step-by-step process is broken, since we obtain a condition on 
the coefficient (1 = which we should be able to consider a parameter here. 

At this point, we have proven that integral distributions within the Vessiot distribution 
exist if, and only if, the algebraic conditions ([33]) are solvable, and that this is equivalent to the 
augmented rank condition (|4ip being satisfied. This in turn is the case precisely if the differential 
equation is involutive. Thus we have characterised the existence of Vessiot connections for 
Vessiot's [28] step-by-step approach. 

8 The existence theorem for flat Vessiot connections 

There remains to analyse the solvability, if we add the differential system (|34p . Its solvability 
is equivalent to the existence of flat Vessiot connections in that each flat Vessiot connection 
of TZi corresponds to a solution of the combined system ([33]) , ([M]) . We first note that the set of 
differential conditions (I34p alone is again an involutive system. 

Proposition 7. The differential conditions ([M]) alone represent an involutive differential equa- 
tion of first order. 

The proof, which is somewhat technical, is given in Appendix[Bj It is based on an evaluation 
of the Jacobi identity for the Lie brackets of the fields Ui. 

If the original equation TZi is analytic, then the quasi-linear system (|34p is analytic, too. Thus 
we may apply the Cartan-Kahler theorem (Theorem [1]) to it which guarantees the existence of 
solutions. The problem is that the combined system ([33]) . ([M]) is not necessarily involutive, as 
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the prolongation of the algebraic equations (j33p may lead to additional differential equations. 
Instead of analysing the effect of these integrability conditions, we proceed as follows. If we as- 
sume that TZi is involutive, then we know from Theorem [2] that the algebraic equations (j33p are 
solvable. Actually, the proof of the Theorem [2] produces an explicit row echelon form of the sys- 
tem matrix given on the right hand side of equation (j4ip . Now we use the interrelations between 

the unknowns ci"''*\ where {a,h) B, to eliminate in the differential conditions (j34p some of 
them by expressing them as linear combinations of the remaining ones. Let I < i < j < n; then 
for all (a, i) where + 1 < a < /3[^^ we find 

j m 

^]°'^ = E E ^7('^iK^'^ + Cf^('^i)' (44) 

while for all {a,i) where + 1 < a < m the correspondence (j35p . given in Remark O holds 
again (see [6l Corollary 3.3.23] for the detailed calculation). In this way, we plug the algebraic 
conditions into the differential conditions and thus get rid of them. 

Example 9. For Example [H where i = 1 and j = 2, we deduce for 

C/i = Xi + Cr + and U2 = X2 + (2^ 3^^ + 8^^ 

from the algebraic condition (137p that —(,2" = ~Ci^ and — = — Ci"^- This is equation (j44p for 
(a, i) = Vx and (a, i) = Wx- The rank conditions ([38]) and ([39|) are satisfied because Vx — wt = 0, 
so that the matrix (^"1 S2 — 0i2) has a vanishing first row. (There is no contraction of 
matrices for j = 2.) There are no interrelations like those in equation p5p because there is 
no index value a for which both i^{du<^) and t*(9ug) would be symbol fields. The differential 
conditions ()34p are 

u^{C2^) - c/2(o = 0, c/i(cr) - f^2(cr) = 0, 

and by substituting ^2"" = ^{'"^ and Ql^"^ = Qi'^ , we can drop the algebraic conditions from the 
system. Since n = 2, there is no further procedure. 

We can now prove the following existence theorem for flat Vessiot connections. 

Theorem 3. Assume that 5-regular coordinates have been chosen for the local representation of 
the analytic differential equation TZi. Then the combined system ()33p . (I34|) is solvable. 

The idea of the proof is this: following the strategy we have just outlined, we eliminate some 
of the unknowns Q. Because of the simple structure of (p4l) . it turns out that we must take 
a closer look only at those equations where the leading derivative is of one of the unknowns we 
eliminate. A somewhat lengthy but straightforward computation shows that these equations 
actually vanish. The remaining equations still form an involutive system. Thus we eventually 
arrive at an analytic involutive differential equation for the coefficient functions Q which is 
solvable according to the Cartan-Kahler theorem. The details are worked out in Appendix [Cl 

Example 10. Consider the first-order equation 

^ . ( Ut = Vt = Wt = Us = 0, Vs = 2Ux + 'iUy, 

^ ' \ = -Ux - 3Uy, Uz =Vx+ 2Wx + 3fy + 4Wy 

in the five independent variables x, y, z, s, t and the three dependent variables u, v, w. It is 
formally integrable, and its symbol is involutive with dimA/i = 8. Thus TZi is an involutive 
equation. For the matrices Si, all of which are 3 x 8-matrices, we find 

/-I ooooo\ /ooo-i 0000 

='1= 0-1000000 , S'o = 0000 -1000 

Vo 0-100000/ Voooo 0-100 
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-0 -1 -2 -3 -4 0\ „ /OOOOOOOOX „ 

-o = 1 -1 0, :r:4 = -2 -4 , ^'i = O'^xS- 

Vo O-l/' * V 1 3 0/' ^ ■^''^ 

For the first two steps in the construction of the fields Ui, the rank conditions are trivially 
satisfied even for the non-contracted matrices. But not so in the third step where we have in the 
row echelon form of the arising 9 x 32-matrix in the 7th row zero entries throughout except in the 
12th column (where we have —2) and in the 17th column (where we have 2). As a consequence, 
we obtain the equality = C2 the rank condition for this step does not hold. However, since 
both Ux and Uy are parametric derivatives and in our ordering Yi = L^{du^) and Y4 = it:{duy), 
this equality is already taken into account in our reduced ansatz and for the matrices the 
rank condition is satisfied. 

Note that the rank condition is first violated when the rank reaches the symbol dimension 
(which is 8). From then on, the rank of the left matrix in (j4ip stagnates at dimTVi while the 
rank of the augmented matrix rises further. The entries breaking the rank condition differ by 
their sign, while their corresponding coefficients in Lemma[2]are collected into one sum and thus 
vanish. Here it shows that contracting the matrices, as explained in Remark [5l is necessary. 



9 Conclusions 

Vessiot's |28j original motivation for the introduction of his theory was to provide an alternative 
proof of the Cartan-Kahler theorem (the same holds for Stormark's presentation [24] of the 
theory). If one takes this point of view, then one may say that our proof of Theorem [3] is 
a "cheat", as it uses the Cartan-Kahler theorem instead of proving it. However, in our opinion, 
this point of view is the main reason why the corresponding proofs in [241 [28] are so difficult to 
read (and actually incomplete as they neglect that the involutivity of the system is a necessary 
condition) . The direct proof of the Cartan-Kahler theorem (for differential equations) given in 
textbooks on the formal theory like |19[I23] is much simpler and more transparent; in particular, 
it makes clear where involution (as opposed to mere formal integrability) is needed. Thus, if the 
only goal consisted of proving such an existence and uniqueness theorem, then there would be 
no need to bother with Vessiot's theory. 

We believe that the real value of Vessiot's theory lies in the fact that it provides via the 
distribution V[7^q] an additional geometric structure on an involutive differential equation TZq 
which is very useful for the further analysis of the equation, that is, after its solvability has been 
established. One possible application is the investigation of certain forms of singular behaviour 
of solutions. In fact, the classical works of Arnold and collaborators (see [2] for an elementary 
introduction) on implicit ordinary differential equations are based on the Vessiot distribution 
(without using this terminology) ; for some further works also in the context of partial differential 
equations or numerical analysis see, for example, |16 [ I25 [ [26]. The basic idea here is that such 
behaviour mainly stems from the fact that at some points the considered involutive distributions 
cease to be transversal to the fibration TZq ^ X so that integral manifolds can no longer be 
interpreted as prolonged solutions in the classical sense. However, it often makes sense to 
consider them a generalised form of (potentially multi-valued) solutions. 

Another application concerns symmetry theory. In classical symmetry theory one always 
considers the differential equation TZq as a submanifold of J^vr and looks then for diffeomorphisms 
of JqTT which are (i) compatible with the contact structure of Jqvr and (ii) leave TZq invariant. 
It was only fairly late realized [1] that this approach yields only what is now called external 
symmetries and that in general further useful symmetries may exist. Using Vessiot's theory, we 
may represent a differential equation as the pair {TZq, V[TZq]), that is, as a manifold together with 
a distribution on it - without any recourse to an ambient space (all relevant properties of Jgvr are 
captured in the Vessiot distribution V[7?.q]). This point of view yields automatically all internal 
symmetries (and is equivalent to the symmetry theory of exterior differential systems). 
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We mentioned already in the Introduction that Vessiot's theory takes an intermediate po- 
sition between the formal theory and the theory of exterior differential systems. It allows for 
the transfer of ideas from the latter to the former one without the need to rewrite a differential 
equation as an exterior system. For example, [231 Section 9.5] contains a formulation of prolon- 
gation structures without exterior forms using the Vessiot distribution insteacH. It should also be 
possible to give an explicit proof of the equivalence of the formal theory and the Cartan-Kahler 
theory on the basis of the results presented here. 

Finally, we comment on the practical side of the Vessiot theory. Throughout this article, we 
made a number of assumptions on the treated differential equation: we assumed that it is first 
rewritten as a first-order system, then that all present algebraic equations are explicitly solved 
and finally that each equation is solved for its principal derivative. While, from a theoretical 
point of view, these operations are always possible under fairly mild regularity assumptions, 
usually they cannot be performed effectively. However, the made assumptions are only needed 
in order to be able to prove results like Lemma [2] on the explicit form of the obstructions to 
involution. These expressions are already bad enough with our simplifying assumptions; for 
more general systems they would hardly be manageable. 

A concrete differential equation need not be transformed to such a special normal form in 
order to apply the Vessiot theory. Even for fully implicit equations of arbitrary order, the 
determination of all integral distributions requires only linear algebra and is easily implemented 
in a computer algebra system. In fact, already Fackerell [5] reported of an implementation of 
Vessiot's theory. Within a recent project thesis [8], Globke implemented large parts of the theory 
in MuPAD. The problems and limitations are here exactly the same as in implementations of 
the Cartan-Kahler theory (see, for example, |10j). 



A Proof of Theorem [2] 

In this section of the appendix, we prove Theorem [2l the existence theorem for integral distribu- 
tions. It is in principle by straightforward matrix calculation and involves a tedious distinction 
of several cases and subcases: we have to compare the entries in the matrix on the right hand 
side of equation ()42p for step j after turning it into row echelon form with the integrability 
conditions and the obstructions to involution as they are given in Lemma [2j For this purpose, 
we fix 1 < i < j and consider the block Si and the entries to its right in the complete matrix. 

Since for a differential equation with an involutive symbol the obstructions to involution va- 
nish and for an involutive differential equation the integrability conditions vanish, too, it follows 
that the augmented rank condition, stated in equation (fl2|) . is equivalent to the differential 
equation being involutive, which in turn is the case if, and only if, the algebraic conditions (j33p 
are satisfied, which is necessary and sufficient for the existence of integral distributions within 
the Vessiot distribution. 

If the order of the differential equation is g > 1, transform it into an equivalent first-order 
equation; for the details of this procedure, see [H Subsection 2.5.1] or [23]. If TZg is involutive, 
then so is TZi (see |23j for a straightforward proof of this) . We assume this first-order equation is 
represented in reduced Cartan normal form (jlOp . We first prove the rank condition (j4ip for the 
homogeneous system. (The proof for the augmented rank condition ()42p follows.) We proceed 
in two steps: We consider the complete matrix at the jth step, 



^Both symmetry theory and prolongation structures were also discussed by Fackerell \5\ via Vessiot theory. 
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(It differs from the matrix on the right hand side of equation (j4ip in that it is not yet contracted.) 
The complete matrix at the jth. step is built from (j — l)j blocks: the stack of j — 1 matrices Si, 
1 < * < J — I, on the left, with each of the Si having another j — 1 blocks to its right, the ith 
of which being Sj and all the others being zero. 

For easier reference, let, for 1 < i, A; < j — 1, be [i, k] the kth. block right of a S^. Then, for 
all 1 < i,k < j — 1, we have 



[z,k] 



" j 5 i — k, 

Omxrj i k. 



For convenience, we set [i,0] := Si. Let, for 1 < g < h < n, 

m^-'' (46) 

denote the matrix that results from writing the block matrices [Si]^ , [Si]^~^^, . . . , [Si]^ (from 
left to right) next to each other. Let [i,k]^ be the hth. block of columns in [i,k], and let, for 
^ 1^ g 1^ h < n, in analogy to the shorthand (j46p . 

[i,kr-^ 

denote the matrix that results from writing the matrices [i,kY, [i,k]^~^^, . .., [i,k]^ (from left 
to right) next to each other. 

If M is any a x 6-matrix, then let ^[M] be the matrix made from the rows with indices (that 
is, labels) c to d, [M]{ the matrix made from the columns with indices e to f and fl^lf 
matrix made from the entries in the rows with indices c to d and in the columns with indices 
from e to /. If the columns of M are grouped into blocks and g denotes which blocks are meant, 
then we write ^[M]^ (with g up right) to show that the first c upper rows are being selected, 
and we write ^[M]^ (with g below right) to show that the last b — d + 1 lower rows are being 
selected. For the block matrix made from M by selecting the rows indexed c to d within the 
block of columns labelled g, we write '^[M]g. This notation is redundant in that the position 
of g does not give new information, but in the calculations to come it increases readability. 

We first turn the complete matrix into row echelon form. Then we contract columns and 
consider the effect. 

The matrices Si are shown in equation pip and have a block-structure of the form 
OsxR — 15 dst 



--: (gij : l<i<A + S,l<j <R + S + T), 



where A, R, S, T are adequate natural numbers or zero, and the blocks are 

{uar ■■ I < a < A,l <r < R) , {has ■■l<a<A,l<s<S), 
{cat ■■ I < a < A,l <t <T) , {dst ■■ I < s < S,l <t <T) , 

the S X S unit matrix I5 or the zero matrix OsxR- (Here the index r is just some index and 
not supposed to be the dimension of any symbol.) Then the substitution 

Cat) ^ ( Cat) + {has) ■ {OsxR "l^ 4i) (47) 

transforms the matrix {gij) into 

/ Qar OAxS Cat + Z]f=l basdst 
\S>SxR -Is dst 
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For the transformation of the complete matrix into row echelon form, we use this obvious method. 
Let 1 < J < n be given. Choose 1 < i < j- Now we eliminate all the non-trivial entries in the 
rows 

f[S,[i,l][i,2]...[i,j-l]] . 

We turn to the block -|^^ [Si] first. Non-trivial entries therein may appear according to equa- 
tion (|3ip only for 1 < h < i in the blocks [Si] . So we fix 1 < h < i and consider the 



block [Si]^. Now set 







(bas) 




(Cat) 


= f [[=,]'.« 


(dst) 




-Is- 


= — 1 (h) , 


OsxR 





[i,l][i,2]...[i,j-l]], (48a) 
_,„^[h,l][h,2]...[h,j-l]], (48b) 



Then it follows that the substitution (1470 leaves {aar) unchanged, turns (bas) into zero as required 
and makes (cat) into {cat + l^f=i ^asC^si)- According to (|48bp the entries X]f=i ^os^^st here 
are of two types: from [S'/i]^,]^ the left part of {dgt), we have the entries in (bag) ■ 

™(h)_^J='/i];j+;^ „, and from ™{h)^J[/i, 2] . . . [h,j - 1]], the right part of (4i), we have those 

of (bas) ■ "(h)^-, [h, I] where 1 < / < j - 1. 

Pl + J- 

Consider the Cat + '^g^i basdst where the factors dst are from the left part of [dst)'- according 
to equation J^ft]/,+i...„ = 0, so 

{bas) ■ ™(h)+J"/i]h+l...n = 0- 

It follows Cat + Yl^s=i basdst = Cat, SO that all entries in [Si]^~^^"'^ ^ the left part of (cgt), remain 

unchanged through the elimination of {b^s) = [Si] . 

For the rest of the proof, we consider the remaining entries Cat + X]s=i basdst , where the 
factors dst are from "(h) [[/i, 2] . . . [h,j — 1]] , the right part of {dgt)- Fix 1 < / < j — 1 and 

consider {has) ■ '^(h)^^[h,l]. 

There are two cases: h = I and h ^ I. For h I we have [h,l] =0 according to the definition 



of the complete matrix ([15|) . thus '"(h) [/i, /] =0 and so 

Pl +1 



{bas) ■'l^(h)^^[hj] = 0. 
Again it follows Cat + X]f=i basdst = Cat, so that for h ^ I all entries in 

f[[i,l][i,2]...[i,j-l]], 
the right part of {cat), remain unchanged through the elimination of {bas) = [^i] , too. 
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For h = I, we have [h,h] = Sj. The structure of Sj, which is (j3ip with i replaced by j, 

Lpl 

fact 



implies that non-vanishing entries are possible in the blocks "(h) [h,h]f^ where 1 < k < j. In 



Jw+^[h,h]k = for l<k<j-l and (49a) 



m 



The entries hasdst for those dgt within \h, h]^ are 



= -1^0). (49b) 



3« 

This matrix has A = rows and T = columns. We consider its entry in row a and 
column t. Setting 7 := + s, 6 := (3^^ + 1 and using S = a[^^ in ([50]) . this entry is 

, (h) 

5^ c!;mcHcP]) (51) 

for all 1 < k < J- Some of the /?f ^ + a^'^^ = m summands vanish because of the special 

entries (09]). As a consequence, from + 1 on, of all the summands C;^(</)")C|((/)J) in ([5T|) at 

(v) 

most one remains: none for k ^ j, exactly one for k = j, namely the one for -f = PI'' +t = 5. 

(i) 

Using the Kronecker-delta, for all /?) + 1 < 5 < m we have 
Now dsn becomes for all 1 < q < /?|*'' and all /?f + 1 < 5 < m 



These are the terms ^5=1 basdst in the case h = I when I < k < j. Now for the subcase 
j + 1 < k < n: these blocks "(h) [h, h]f^ = "^^h) ["j]/;) again on account of the structure of Sj, 

which is (j3ip with i replaced by j, are zero (if they exist at all; they do not for k with a^*'^ = 0). 
So in this case we have 



(bas) 



l...n 



which contains the terms X]j,=i basdst = 0. As a consequence, here, in the case h = I, the 
terms Cat + X^f=i basdst are of the following form. The Cat of interest in line (j48ap are those 

in the block [^,^]. For h i, we have [i,h] = 0, thus [i, /i] = 0. Non-trivial entries Cat 

are possible only for h = i; m. that case = Sj, thus [i, /i] = ["j]- According to the 

structure of Sj, which is (j3ip with i replaced by j, the only entries Cat within the block [Sj\ 
which may not vanish are, for 1 < k < j , those of the form — C|((/>"). They make up the block 
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k 





1 


2 


i - 1 


i 


i + 1 ... j -1 


j 




1 


d.[d2) 














2 




3.(d2) 


3.(dl) 


3.(c) 


3.(b) 


3.(a) 


h 




















i - 1 






3.(d2) 










i 








2.(a) 


2.(b) 



Figure 2. Possible combinations of h and fc in the terms u^j. used as labels for columns in the contracted 
matrix. Each term w^j. corresponds to one subcase in the consideration of cases 2. h = i and 3. h < i. 
Each of these subcases refers to one of the terms in square brackets in Lemma [2] as follows: 2. (a) - 
line (ini), 2.(b) ~ line HID, 3. (a) - line 3.(b) - line (HH), 3.(c) - line 3.(dl) - line 3.(d2) - 
line (fT3l) . 



So for any row index a and any column index 6 = j3^^ + i we have Cat = ~^hiCs{4''j) as the 
most general form of an entry. 

Since (|5Hp is Yls=i basdst, it follows that Cat + Y,s=i basdgt is 



-5f,,cHc^J)+ J2 C!;{cP2)CHcP]) + 6,jCnct>f). (52) 

(i) (k) 

bmce 1 < a < PI'' and Pf' + 1 < 5 < m, any term in the row echelon form of the complete 
matrix without a —1 in a negative unit block somewhere to its left has this form or is a zero 

in [Sif for some i < j (in which case it remains zero throughout the elementary row 
transformations and so does not influence the rank). This means, if all these expressions vanish 
(when contracted), the rank condition is satisfied. To show that they do vanish (when con- 

tracted) for a system with an involutive symbol, we consider them as the new entries in [i, h] 
with 1 < h < j . Now with regard to the relation between h and i there are three cases: h > i, 
h = i and h < i. We consider them in that order. 

1. Let h>i. Then according to the structure of Si, all the C!^{(pf) and C^(0") in (f52D 
vanish. Since h ^ i, 6hi = 0, thus all of ([52]) vanishes. 

2. Let h = i. We shall consider several subcases for h = i, and for h < i after that, which 
may be labelled by second-order derivatives where h and k are combined as shown in 
Fig. [2j For fixed 6, any belongs, according to its indices h and k, to exactly one of the 
blocks marked by the case denominations 2. (a), 2.(b) and 3. (a) to 3.(d2). 

It turns out that (|52p is the common form of all the sums that appear in the squared 
brackets of Lemma [2] as the coefficients of second-order derivatives n^^. Not only can the 
case distinctions of the following argument, 2. (a), 2.(b) and 3. (a) to 3.(d2), be labelled by 
these u^f^f., but in fact a case labelled n^^ is dealt with by using the fact that according 
to Lemma [2] the coefficient in square brackets of that same vanishes for an involutive 
system. The correspondence of the cases and the terms u^^, is given in the caption of 

Fig.m 

The cases h = i with subcase k < i and h < i with subcase k < i need not be considered 
because the columns of the complete matrix are to be contracted when proving the rank 
condition. This contraction concerns those columns of the complete matrix which are 
labelled by the same second-order derivative, and each second-order derivative is used 
exactly once in the argument. 
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(a) For the first subcase of /i = i let i < A; < j. Then ([52]) becomes 



(k) 

where /3| + 1 < 5 < m. This vanishes for a system with an involutive symbol according 
to Lemma [21 line (jl7|) . 

(b) For the second subcase of h = i choose k = j. Then (|52p becomes 



where + 1 < 6 < m. According to Lemma [21 line p^ . this vanishes for a system with 
an involutive symbol. 

3. Let h < i. We consider several subcases with regard to the relation between h < i, j and k. 
(a) First choose k = j. Then (|52p becomes 



5^ C,^<^nC^j(<^J) + C7,^ 



where Pi'^ + I < 6 < m. According to Lemma [2l line (jlSp . this vanishes for a system with 
an involutive symbol. 

(b) For the second subcase choose i < k < j. Then ()52p becomes 



7=/3f'+l 

where /jj'^^ + 1 < 5 < m. According to Lemma [21 line ([16]) . this vanishes for a system with 
an involutive symbol. 

(c) For the third subcase choose k = i. Then (|52p becomes 

E C,^<^?)C7K<A]) (53) 

7=/3l'''+l 

where PI ' + 1 < 5 < m. Since we have h < k = i, for any such 5 the cross derivative 



u, 



ih 



u^. labels two columns in the complete matrix: the one with label 6 in [i, /i]* and 



the one with label 5 in [h^]'^ which according to ([52p has the new entries 



7=/3f'+l 
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where 15^^^ + \ < 5 < m. In the current subcase, ([5^ becomes 

-Cl{^'^)+ E C;mcl{<P,) (55) 

where /jf ^ + l<5< m. Since the columns of both ([53]) and ([55]) are labelled by the same 
second-order derivatives, we have to contract them, which means adding their new entries. 
For all /jf ^ + 1 < (5 < (3^ this yields 



which, for a system with an involutive symbol, vanishes according to Lemma[2l line (|15ap . 

and for all pf^ + 1 < 6 < m 



which, for a system with an involutive symbol, again vanishes according to Lemma [21 
line (fTSb]) . 

(d) For the fourth subcase choose k < i. Under this assumption, we have to distinguish two 
further subcases: k = h and k < h. The subcase k > h need not be considered since 
'^kh ~ ''4ik second-order derivative is used only once to label a column in the 

contracted matrix. 

(dl) First consider k < i and k = h. Then still h < i, furthermore k < j, and ([52p becomes 
7=/3j'"+l 

where /jj'*^ + 1 < 5 < m. According to Lemma [21 line (|13p . this vanishes for a system with 
an involutive symbol. 

(d2) At last consider h < k < i. Then k < j, and ([52]) becomes 

^ C,^<^r)C|('^J) (56) 

7=/3l'''+l 
(k) 

where /?)^ + 1 < 5 < m. Since we have h < k < i, ior any such 6 the cross derivative 
^fc/i = Uf^ij. labels two columns in the complete matrix: the one with label 6 in [z, /i] 
and the one with label 5 in -^^ [i, A;] which according to (|52p has the new entries 

-6,,C^i^-)+ C!;{c^?)C^{4>]) + 6h,Cl{^t) (57) 
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where [5^^^ + 1 < 5 < m. In the current subcase, ([57|) becomes 

c!;{4>t)C^{<p]) (58) 

where /?f ^ + l<5< m. Since the columns of both (j56p and (j58p are labelled by the same 
second-order derivatives, we have to contract them, which means adding their new entries. 
This yields 



for Z?!'^'* + 1 < (5 < /?|'^'' and vanishes according to Lemma [21 line (I14ap . and 

for Z?!'^^ + 1 < 5 < m and vanishes according to Lemma [21 line ()14bp . 

What we have shown is that in the contracted matrix all the entries without some entry — 1 
in a negative unit block to their left may be eliminated by elementary row transformations if 
the equation has an involutive symbol. Thus under this assumption the rank condition (j4ip is 
satisfied. 

Now for the augmented rank condition, equation ()42p . To transform the augmented complete 
matrix into row echelon form, we use for each j and i the same procedure as for the transfor- 
mation of the non-augmented complete matrix, except that now the matrices (cat) and (dst) 
are augmented by one more column each as follows. Fix I < j < n. Let 1 < i < j. Then, 
according to the structure of the augmented complete matrix, for its transformation into row 
echelon form we have to eliminate for 1 < h < i the entries in the blocks [S'j]'^, as we did for the 



non-augmented complete matrix given in equation (j45p . We have to consider the effect of these 
transformations on the additional entries which make up the rightmost column in the augmented 
complete matrix. These are —G>jj, • • • 5 given in equation ([28l) . Of these entries, only 

8^^^ (i) 

—G>jj, • • • 5 affected (since we eliminate the entries which are in rows 1 to 

of the matrices [S'j]'^). We add them as the rightmost column in the augmented matrix (cat)- 
Now fix 1 < h < i. Then augment the matrix {dgt), used in the process of eliminating the 

entries in [Si] , by adding the entries ' ~^hj ' • • • > ~^h!j ^ rightmost column, 

in accordance with the structure of the augmented complete matrix. Then the substitution (j47p 
yields as the transformed entries Cat + Sf=i basdst the same as for the transformation of the 
non-augmented matrix except, of course, for the new last column. Now let denote t the index 
of this last column. Fix some row index 1 < a < /3{ . Then the entry Cat = ~^fj transforms as 
follows: For all 1 < /i < i we have to add to it 
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Since here the matrix {h^'y) = [^i]^ , this is the product of the row with index a in the matrix [S 

and the transpose of the vector i~^hj i ■ ■ ■ i~^h^j)- According to the structure 

of [^i]^ as defined in equation ([5D|) . this equals 

m 

The entries 0]^ - can be taken from equation (p8|) . Since < 7 < m, we have O^- = C^^^((/)^) 

if /?f ^ + 1 < 7 < /^P and Ol - = if /J^ + 1 < 7 < m. Therefore we get 



(i) 

as the summand for each h to be added to Cat = —&fi- Since we have —0Z 

—C\^\(j)J) + C^p {4)f). Therefore the entry Cat transforms into 



which is the integrabihty condition in hne p2p of Lemma [21 except for the sign. The difference 
in sign comes from the fact that the augmented complete matrix describes the system of equa- 
tions (j40|) . where the entries from the row with index a in the matrix Si are on the opposite side 
from the entries from the row with index a in the matrix Sj and the inhomogeneous term —O'iji 
while in Lemma [2] for each i, j and a the corresponding equation is set to zero, if the differential 
equation is involutive. 

This means that the augmented rank condition holds for all 1 < i < j < n if, and only if, the 
equation has an involutive symbol and is formally integrable. 



B Proof of Proposition [7] 

We are going to show that the system (j34p is involutive. Consider the basis {Ui : 1 < i < n) 
of the distribution U given by Ui = Xi -\- Cf ^ with yet undetermined coefficient functions 
Cf S J-{TZi). For all 1 < i < n and 1 < p < r, the independent variables of the functions 
are the coordinates on TZi, which are x, u and all such that {a,h) B. To apply a vector 
field Uj = d^j + • • • to a function (^f includes a derivation with respect to . We order the 
independent variables such that if j > i, then is greater than x*, and each is greater than 
all the variables and where (a, h) B. For any equation Hfj within the system (p4]l . the 
application of the vector field Uj = d^j + • • • to yields d(f /dx^ as the leader of that equation; 
therefore equation Hfj is of class j, and the equations of maximal class are Hf^; the equations 
of second highest class in the system are H^n-i ^'^'^ There are only equations of a 

class indicated by some index 2 < j < n. 

From the Jacobi identity for vector fields Ui, Uj and Uk where 1 < i < j < A; < n, we have 

[U„ [Uj,Uk]] + [Uj, [Uk, [/.]] + [Uk, [Ui, Uj]] = 0. 
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The structure equations ([32]) for the vector fields Uh and the definitions of G^j and Hfj in 
equations (|33|) and (|34|) imply that this is 



= G%{U„ ZJ + U,{G'j,)Z, + H^,[Ui, Yp] + U,{H^,)Yp 

+ GUUj, Z,] + Uj{Gl,)Z, + H^U^Yp] + 

+ Gijpk, ZJ + C/fc(G^,)^c + Hf^[Uk,Yp] + Uk{HP)Yp. 

The combined system (j33l [31]) means that all G^^ = and all H^^ = which implies that U is 
involutive, which it is, being in triangular form, exactly if all [Ua, Ub] = 0. This leaves only 

= mG'j,) + [/,(GL) + UkiG'^^)}Z, + mH^,) + U,{Hl^) + Uk{Hf^)}Yp. 

As part of a basis for V'[7^g], the vector fields Zc and Yp are linearly independent, which means 
their coefficients must vanish individually. So in particular 

Under the assumption i < j < k, the term Uk{Hf-) contains derivations with respect to 
of Ui{(^j) and Uj{Cf). Thus, according to our order, this is a non-multiplicative prolongation, 
and the remaining terms are multiplicative prolongations. But since any non-multiplicative 
prolongation within the system (j34p must be of such a form, it is a linear combination of 
multiplicative prolongations. Therefore, no integrability conditions arise from cross-derivatives 
(and none arise from a prolongation of lower order equations since all equations of the system 
are of first order). 

If we set dC,f /dx^ =: (Cf )j for the leaders and 

and solve each equation of the system ()34p for its leader, then it takes the form 

(Cf)n = C/2(CD-f>n(Cf), 



(cnn-i = c/i(CLi)-c^«-i(Cf), 



(Cf; 


>n-l = u^ic. 


-1) - Un- 


-l(Cf), 
-f>n-l(CL2 






2(cLi) - 


(Cf; 


>3 = U2{0 - 


f>3(Cf), 




(Cf; 


i3 = c/i(Cf)- 






(Cf; 


i2 = c/i(Cf)- 


f>2(Cf); 





MT) 



here for each line I < p < r. Therefore the system ()34rp is in Cartan normal form given in 
Definition [71 

Now one can prove (see O Lemma 2.4.29]) that such a differential equation TZi is involutive 
if, and only if, all non-multiplicative prolongations of the equations (|8ap - ([8c]) and all formal 
derivatives with respect to all the x* of the algebraic equations ([8dp are dependent on the 
equations of the system ([5]) and its multiplicative prolongations only. Therefore, it follows that 
the system ([3l[) is involutive. 
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C Proof of Theorem [3] 

We consider the combined system of algebraic and differential conditions (|33p . (j34|) and want 
to show that it has a solution. We follow the strategy outlined above and eliminate some of 
the unknowns C]*. As we consider each of the equations of (j34p as being solved for its derivative 

dc!f'^'^/dx^ of highest class j, as given in equation (j34rp . we must take a closer look only at those 
equations where this leading derivative is of one of the unknowns we eliminate. The structure of 
the vectors Cij given in equations (j35p and (j44p . shows which ones these are. Let k be such that 
2 < k < n. Then for the subsystem of the equations of class k in the system ([M]) . the equations 
which hold the following terms are concerned: 



here, for any Uk{Q^^''^^), we have /Sj'^^ + 1 < (5 < m. We now show that these equations vanish. 
The proof is by straightforward calculation, though tedious and requiring a case distinction. Let 
1 < i < k. Fix some /jj^^ + 1 < /? < m. Consider the equation 

Ur{C'i^'^)=Uu{Cf''^). (59) 

Then h < i < k. According to the structure of the vector the entries of which in its hth block 
are of two kinds, there are two cases. 

1. The interrelation for (^i^'^^ is an equality: Ci^''^^ = Cjf'^^- This is so if, and only if, 
(i) 

PI + 1 < P < m according to the structure of Q. Now there arise two subcases, 
(a) The other interrelation is an equality, too: cf '''^ = Ch'''\ This is so if, and only if, 

(k) 

Pj^+l<P<m according to the structure of Ck- In this subcase, equation (j59p 
becomes 

U^{C^^) = U,{ct'^). (60) 

Since the system (fM]) contains the equalities Ui[Cff'''^) = Uh{Ci^'''^) ^'^'^ ^fclcl^'*"*) ~ 
Uh{Cjf'^^), equation ([60]) becomes 

Since i < k and ^[''^ + 1 < /? < from the structure of Ck follows Ci^''^^ = Ci^'^^ ■ 



Thus, equation (j59j) vanishes, 
(b) The other interrelation is an affine-linear combination: 

A; m 

(i) (k) 

This is so if, and only if, /3| + 1 < /3 < according to the structure of ■ In this 
subcase, the term Ui{cjf''^^) in equation ([59j) becomes 

k m 

mcI'''^)=y: e (61a) 
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+E E HcM))^^r:^^^+Hc'^'m- (61b) 



The term Uk{c!f'^'^) in equation ([59]) becomes 



E E ^."(^f)f^^(C^"^) (62a) 



+E E Mc^M))cr' + Mcl'\4))-, (62b) 

here we have the first equahty because we are considering the first main case, the 
second equahty because of the structure of the system (|34p and the third equahty 

(i) (k) 

according to the structure of Cfc > since i < k and because /3J + 1 < f3 < f3i . 



Substituting (|6T|) and (|62]) in equation (|59l) and factoring out, we get 

= E E ^^7(^f){f^(d""^)-f^4Ci^"^)} (63a) 



+E E mcM))cr'-McM))crn (63b) 

+ f^^(^^i'^(^f))-^.(C^fH^f))- (63c) 

Line (j63cp contains the Lie bracket [J7j, C^^^]((^^). According to the structure of the 
system (pll) . the term (I63ap vanishes. If the terms (|63bp and ()63cp vanish, too, then 
so does equation (j34p . Otherwise they form a new algebraic condition for (j34p . which 

can be solved for some function Cf!^'^^ ■ Substituting this function in (I34p does not 
change the classes or the numbers of the single equations therein. Thus, equation (j59p 
vanishes. 

2. The interrelation for cP'^^ is an affine-linear combination: 



cP^ = E E q(^f)c^"^ + Gi^H^f) 



This is so if, and only if, + Pi^ according to the structure of Q. Since we 

(i) (k) 

have h < i < k and /3J 1^ Pi , according to the structure of the other interrelation is 
an affine-linear combination, too: 
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Thus, equation ([59|) becomes 



^ = E E c^M)Mc^r^) -Y. E cM)Hct'^) (64a) 



+ E E Mcm)c'r'-i: E HcMM'' (64b) 



+ f^4C^i'H^f))-f^4^^i'^(^f))- (64c) 

In part dMSD, the terms f/fc(Cl'^'''^) and f/»(cf'*^) are equal to UhiCk'''^) and ?7/^(ci^'''^) 
according to the structure of the system (j34|) . Thus, equation (j64p becomes 



k m 



+ (l64bl) + dMcD. 



Factoring out the vector field in part ()64aip . this equals 



= ^ME E ^7(^f)Cr"^-E E C7.^(^SCP1 (65a) 



-EE ^^(c^7(^f))cr^-E E Mctmcr] (65b) 

+ (IMlo) + (IMt). (65c) 
Now one can show (see [6l Corollary 3.3.25 (for j = k)]) that the term (|65ap equals 

f^.(^^f^(^f)-^^I'^(<^f)), 

which does not contain any C^''^'"^ or Qf'^'^ any more; it is an algebraic expression instead of 
the differential expression that it seems to be when written in the form (j65ap . The other 
terms, ()65bp and (|65cp . are algebraic, too. So all of equation has shown to be an 
algebraic condition when the interrelations between the entries of the vectors Qh-, d and Qt-, 
as noted in equations ([35]) and (jM]), are taken into account. 

If this new algebraic condition for the system (|34p vanishes, equation (|59p vanishes. Oth- 
erwise, this new algebraic condition given in equation (|65p now appears as 

^ = Mc\'\4)-c^\4)) S) 

/ i m km \ 



-EE Mc^M))cr'-Y. E McM))cr] mi 

i m k m 

+ E E f^^(^7(^f))cr"^-E E u.ic^MM'^ m 

+Uk{ci:\^^))-u^{ci^\4))- m 
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Collecting terms in lines (|65al ) and (|64cp . this yields 
= pibl) + pbl) 

The lower line contains the Lie brackets [Uh,cl^^]{(j)^) and [Uk,cl^^]{(/)^). There must be 

some non- vanishing summand containing a factor C^'^''^'', Ci^''^\ C^'"^ or Ch''^^ ■ As we did in 
case 1. (b), we solve (|64|) for this non- vanishing factor and substitute it into the system (|34p . 
which does not change the class of any equation therein. Therefore equation (j59p drops 
out from the system (|34p . 

Now we have shown that all those equations vanish where the leading derivative is subject 
to being substituted through the interrelations concerning the coefficient function Cfc^''^ In the 
system (j34rp . these are the equations with the leaders 



V'fc-l /fc' \^k-l Jk' \^k-l Jk' 

here 2 < k < n and p[^'^ + 1 < P < m. The remaining equations still form an involutive system 
(we may enumerate the remaining (^f in such a way that no gaps appear) as the considerations 
for the system ()34p in Proposition [7] apply likewise. Thus we eventually arrive at an analytic 
involutive differential equation for the coefficient functions which is solvable according to the 
Cartan-Kahler theorem (Theorem [1]) . 
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